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13  ABSTRACT  (Maximum 200 words)  A  momentum  exchange  controller  where  the  time  rate  of  change 
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levered  flexible  appendage  undergoing  a  single  axis  maneuver;  and  (2)  as  a  component 
of  an  adaptive  feedback  control  strategy  for  the  retrieval  of  an  orbiting  tethered 
antenna/reflector  system.  The  optimal  linear  quadratic  Gaussian  (LQG)  digital  con¬ 
trol  of  the  orbiting  tethered  antenna/reflector  system  is  analyzed.  The  flexibility 
of  both  the  antenna  and  the  tether  are  included  in  this  high  order  system  model. 

With  eight  point  actuators  optimally  positioned  together  with  tether  tension  modula¬ 
tion  it  is  seen  that  the  degree  of  controllability  is  very  high.  A  method  of  mea¬ 
suring  tether  transverse  motions  is  proposed  and  is  required  to  guarantee  system 
observability.  An  analytical  formulation  for  the  modelling  of  the  solar  radiation 
disturbance  on  the  tethered  antenna/reflector  system  is  introduced.  A  control  law 
based  on  tension  modulation  where  the  gains  are  based  on  the  linear  quadratic  regu¬ 
lator  theory  is  able  to  maintain  satisfactory  pointing  accuracy  for  low  and  moderate 
altitude  orbits  under  the  influence  of  solar  pressure.  For  the  higher  altitudes  a 
combination  of  tether  tension  modulation  together  with  active (actuator)  control  will 
be  required. 
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Summary 


A  momentum  exchange  controller  where  the  time  rate  of  change 
of  the  flexible  momentum  relative  to  the  rigid  body  motion  is  used 
as  a  part  of  a  feedback  control  law  for  maneuvers  and  flexible 
vibration  suppression  is  introduced.  This  control  concept  is 
applied  to:  (1)  a  model  of  a  rigid  hub  (base)  with  a  cantilevered 
flexible  appendage  undergoing  a  single  axis  maneuver;  and  (2)  as  a 
component  of  an  adaptive  feedback  control  strategy  for  the 
retrieval  of  an  orbiting  tethered  antenna/reflector  system.  The 
optimal  linear  quadratic  Gaussian  (LQG)  digital  control  of  the 
orbiting  tethered  antenna /ref lector  system  is  analyzed.  The 
flexibility  of  both  the  antenna  and  the  tether  are  included  in  this 
high  order  system  model.  With  eight  point  actuators  optimally 
positioned  together  with  tether  tension  modulation  it  is  seen  that 
the  degree  of  controllability  is  very  high.  A  method  of  measuring 
tether  transverse  motions  is  proposed  and  is  required  to  guarantee 
system  observability.  An  analytical  formulation  for  the  modelling 
of  the  solar  radiation  disturbance  on  the  tethered 
antenna/ reflector  system  is  introduced.  A  control  law  based  on 
tension  modulation  where  the  gains  are  based  on  the  linear 
quadratic  regulator  theory  is  able  to  maintain  satisfactory 
pointing  accuracy  for  low  and  moderate  altitude  orbits  under  the 
influence  of  solar  pressure.  For  the  higher  altitudes  a 
combination  of  tether  tension  modulation  together  with  active 
(actuator)  control  will  be  required. 


ACKNOWLEDGEMENT 


This  research  has  been  conducted  under  the  direction  of  Dr. 
Spencer  Wu,  Program  Manager,  Aerospace  Sciences,  Air  Force  Office 
of  Scientific  Research,  Bolling  AFB,  Washington,  D.C.  Appreciation 
is  expressed  for  the  strong  encouragement  and  useful  comments 
provided  by  Dr.  Wu  during  the  course  of  this  study.  Thanks  are 
also  extended  to  Dr.  Wu  for  making  time  available  for  an  oral 
progress  report  by  all  participants  during  his  visit  to  Howard 
University. 


TABLE  OF  CONTENTS 


Stunmary  ill 

ACKNOWLEDGEMENT  i  V 

1 .  INTRODUCTION  1 . 1 

1.1  Feasibility  of  Concept  based  on 

Existing  Work  1.1 

1 . 2  Relevance  to  SDI  and  Military  1 . 3 

Multibeaia  Communications 

1.3  Outline  of  the  Research  Reported  1.3 

2.  MOMENTUM  EXCHANGE:  FEEDBACK  CONTROL  OF  FLEXIBLE 

SPACECRAFT  MANEUVERS  AND  VIBRATION  2.1 

2 . 1  Introduction  2 . 1 

2 . 2  Momentvun  Exchange  Feedback  Control  2 . 3 

2.3  Independent  Flexible  Control  for 

Elastic  Vibration  2 . 7 

2 . 4  Large  Angle  Maneuvers  2 . 8 

2 . 5  Conclusion  2 . 9 

3.  ADAPTIVE  FEEDBACK  CONTROL  FOR  THE  RETRIEVAL 

OF  AN  ORBITING  TETHERED  ANTENNA/REFLECTOR  SYSTEM  3.1 

3 . 1  Introduction  3 . 1 

3 . 2  The  System  Equations  of  Motion  3 . 2 

3 . 3  Momentum  Exchange  Feedback  Control  for  System 
In-Plane  Motion  at  a  Fixed  Tether  Reference 

Length  3 . 8 

3 . 4  Retrieval  3 . 9 

3.5  Numerical  Results  3.11 

3.6  Conclusion  3.12 

4.  THE  OPTIMAL  LQG  DIGITAL  CONTROL  OF  AN  ORBITING 

ANTENNA/REFLECTOR  SYSTEM  4.1 


4 . 1  Introduction  4 . 2 

4.2  Mathematical  Model  4.3 

4 . 3  Simulations  4 . 14 

4 . 4  Conclusions  4 . 16 

5.  EFFECTS  OF  SOLAR  RADIATION  PRESSURE  ON  THE 

TETHERED  ANTENNA/ REFLECTOR  SUBSATELLITE  SYSTEM  5.1 

5 . 1  Introduction  5 . 1 

5.2  Modelling  of  the  Solar  Radiation 

Disturbance  Torques  5 . 4 

5.3  Modeling  of  the  System  Equations  of  Motion  5.12 

5.4  Numerical  Simulations  5.16 

5.5  Conclusions  5.17 

6.  CONCLUSIONS  AND  RECOMMENDATIONS  6.1 


1. 


INTRODUCTION 


1.1  Feasibility  of  Concept  based  on  Existing  Work 

In  an  earlier  paper  [1]*  Yasaka  suggested  that  a  large  Earth- 
oriented  flexible  orbiting  antenna  system  culd  be  stabilized  by 
using  the  gravity-gradient  torques  through  the  adjustment  of  the 
lengths  of  a  series  of  tethers  which  would  be  connected  between  the 
feed  panels  and  the  reflector  surface.  The  structural  feasibility 
and  a  related  stress  analysis  was  presented.  In  a  subsequent 
related  treatment,  Bainum  and  Kumar  [2]  investigated  the  dynamics 
and  stability  of  an  orbiting  flexible  shallow  spherical  shell  with 
a  dumbbell  connected  to  the  shell  at  its  apex  by  a  spring  loaded 
double  gimballed  joint.  In  addition  to  providing  a  favorable 
composite  moment  of  inertia  distribution  for  gravitational 
stabilization,  through  the  use  of  internal  magnetics,  (eddy 
current)  passive  damping  can  also  be  utilized.  Subsequently, 
Bainum  and  Reddy  [3]  have  considered  the  shape  and  orientation 
control  of  the  shell  antenna  by  including  active  control  elements 
in  addition  to  the  passive  damping.  Numerical  results  verified  the 
significant  savings  in  fuel  consumption  that  could  be  realized  when 
the  hybrid  she 11 -dumbbell -actuator  system  is  acting  together  as 
compared  with  the  active  thrusters  operating  alone. 

In  the  first  year's  contractual  effort  [4]  the  basic  equations 
of  motion  and  related  stability  analyses  for  a  tethered  antenna 


‘For  references  sited  in  this  report  please  refer  to  reference 
list  at  the  end  of  ech  chapter. 


reflector  system  during  station  keeping  were  developed,  and  optimal 
tension  control  laws  for  inplane  motion  control  were  synthesized. 
It  was  concluded  that  additional  types  of  control  would  be  required 
for  out-of-plane  motion  damping  and  also  for  providing  shape 
control  of  the  flexible  antenna  surface. 

A  comprehensive  review  of  the  literature  related  to  the 
control  of  various  classes  of  tethered  satellite  systems  for 
deployment,  station  keeping,  and  retrieval  was  conducted.  Both 
linear  and  nonlinear  control  laws  were  reviewed,  with  the 
conclusion  that  control  laws  based  on  optimal  control  theory  offer 
the  greatest  potential  for  the  proposed  orbiting  tethered 
ref lector/antenna  systems. [4] 

The  second  year's  research  [5]  concentrated  on  the  use  of 
active  control  thrusters  to  provide  first  order  control  for  the 
roll/yaw  motions  of  the  shell  and  also  the  out-of-plane  tether 
swing  motion.  Actuators  would  be  positioned  at  the  end  of  the 
rigid  boom,  through  which  the  tether  is  deployed  and  connected  to 
the  shell,  and  also  placed  along  the  periphery  of  the  shell. 
Additional  thrusters  mounted  normal  to  the  shell's  largest  surface 
could  also  be  used  in  combination  with  tether  tension  modulation  to 
ensure  both  orientation  and  antenna  shape  control. 

Control  laws  were  also  developed  that  could  be  used  for 
deployment  of  the  tethered  antenna  system  where  the  actual  control 
gains  could  be  adjusted  in  a  piecewise  adaptive  manner  depending  on 
the  difference  between  the  commanded  and  true  length  and  the  mass 
of  the  deployed  tether.  A  deployment  strategy  was  also  developed 
which  could  be  adopted  to  an  orbiting  test  scale  model  of  the 


system.^**  Finally,  some  preliminary  results  were  presented  'or  the 
defensive  near-minimum  time  slewing  maneuver  of  the 
antenna/reflector  with  the  tether  either  fully  retrieved  or 
(assumed)  separated.  The  two-point  boundary  value  problem 
resulting  from  the  application  of  Pontryagin's  maximum  principle 
was  (numerically)  solved  by  the  quasilinearization  algorithm.**' 

1.2  Relevance  to  SOI  and  Military  Mutibeam  Communications 

Associated  with  the  capability  to  orient  a  large  flexible 

antenna/reflector  type  of  device  accurately  while  at  the  same  time 
maintain  the  surface  shape  to  within  centimeters  or  even 
millimeters  are  many  applications  in  both  the  military  and  civilian 
fields.  For  example,  high  energy  beams  can  be  generated  by  a  power 
source  and  reflected  from  specific  known  points  on  the  reflector 
surface  to  preselected  targets.  In  the  very  important 
communications  field,  such  an  antenna  surface  can  receive  multibeam 
communication  waves  from  electronic  feed  devices  and  transmit  these 
to  a  variety  of  small  mobile  receivers  to  comprise  strategic 
communication  links  during  early,  critical  phases  of  an  attack  when 
larger,  fixed  land-based  antennas  would  be  far  more  vulnerable  to 
observation/damage.  Such  devices  could  also  be  employed  to 
transmit  coded  electronic  mail  rapidly  over  different  communication 
channels. 

1.3  Outline  of  the  Research  Reported 

The  second  chapter  introduces  the  concept  of  a  momentum 
exchange  controller  where  the  time  rate  of  change  of  the  flexible 


momentum  of  une  system  relative  to  the  rigid  body  motion  is  used  as 
a  part  of  the  feedback  control  law  for  maneuvers  and  vibration 
'suppression.  A  particular  example  of  a  rigid  central  part 
connected  to  a  pair  of  cantilevered  flexible  appendages,  such  as 
solar  panels,  is  considered.  Lyapunov  stability  techniques  are 
used  to  obtain  upper  and  lower  bounds  on  the  control  law  gains. 
With  the  presence  of  this  type  of  feedback  control,  an  additional 
independent  control  system  acting  on  the  flexible  parts  can  be 
designed  for  further  vibration  suppression.  A  paper  based  on  this 
research  was  presented  at  the  AAS/AIAA  Astrodynamics  Specialist 
Conference,  Durango,  Co.,  Aug.  19-22,  1991. 

The  next  chapter  also  uses  the  momentum  exchange  control 
technique  as  a  part  of  an  adaptive  feedback  control  concept  for  the 
retrieval  of  the  orbiting  tether  antenna/ ref lector  system.  The 
time  rate  of  change  of  that  part  of  the  angular  momentum  due  to 
the  tether  motion  from  the  local  vertical  is  used  as  a  part  of  the 
feedback  law  for  the  antenna  attitude  control  system.  This  is  used 
in  conjunction  with  an  open-loop  exponentially  decreasing  tether 
commanded  length  control  during  retrieval.  The  selection  of  the 
dimensionless  control  gains  is  b' sed  both  on  system  stability 
considerations  when  the  system  is  linearized  about  a  given 
reference  length,  and  on  numerical  simulation  of  the  nonlinear  time 
varying  system.  The  feedback  control  gains  are  then  adapted  to  the 
instantaneously  changing  tether  length.  This  chapter  is  based  on 
results  presented  in  a  paper  at  the  Second  Joint  Japan/U.S.A. 
Conference  on  Adaptive  Structures,  Nagoya,  Japan,  Nov.  12-14,  1991. 

Chapter  4  considers  the  optimal  LQG  (linear  quadratic 


Gaussian)  digital  control  of  the  orbiting  flexible  tethered  antenna 
system.  The  flexibility  of  both  the  tether  as  well  as  the  antenna 
(shell)  is  included.  A  comprehensive  computer  simulation  program 
is  developed  to  also  include  the  important  effects  of  state 
measurement,  plant  and  measurement  noise,  system  controllability 
and  observability.  The  LQG  technique  is  employed  to  synthesize  a 
controller  with  the  Kalman  filter  to  deal  with  measurement  and 
plant  noise  in  the  presence  of  sampled  data  observations.  The 
optimum  number  and  location  of  the  actuators  are  determined  by 
means  of  the  concept  of  degree  of  controllability  and  related 
numerical  simulations.  The  results  contained  in  Chapter  4  willl 
also  be  included  in  a  paper  accepted  for  presentation  at  the 
AAS/AIAA  Space  Flight  Mechanics  Meeting  in  Colorado  Springs,  Co., 
Feb.  24-26,  1992. 

In  the  following  chapter  the  effect  of  solar  radiation 
pressure  on  the  tethered  antenna/reflector  subsatellite  system  is 
considered.  The  solar  radiation  pressure  is  one  of  the  dominant 
sources  of  disturbances,  especially  in  the  higher  altitude  orbits. 
One  of  the  objectives  of  Chapter  5  is  to  develop  the  analytical 
formulations  for  the  modelling  of  solar  radiation  disturbance 
torques  on  the  tethered  antenna  subsatellite  system.  A  second 
objective  is  to  evaluate  the  effects  of  these  disturbances  (for 
station  keeping  operations)  to  determine  at  what  altitude  levels 
the  tension  control  law  is  still  able  to  maintain  the  required 
pointing  accuracy,  and  to  develop  a  hybrid  compensation  control 
strategy  (involving  tether  tension  plus  some  kind  of  actuators 
attached  to  the  antenna)  for  the  situations  where  the  tension 


modulation  control  alone  can  not  meet  the  mission  pointing 
accuracy.  Papers  based  on  this  chapter  (Chapter  5)  have  been 
offered  for  possible  presentation  at  the  1992  AIAA/AAS 
Astrodynamics  Conference,  Hilton  Head,  S.C.,  Aug.  10-12,  1992  and 
also  at  The  World  Space  Congress  ( lAF/ COSPAR) ,  Washington,  D.C., 
Aug.  28-Sept.  5,  1992. 

Chapter  6  summarizes  some  concluding  remarks  and  suggestions 
for  further  continuing  research. 
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1.6 


2.  MOMENTUM  EXCHANGE:  FEEDBACK  CONTROL  OF  FLEXIBLE 
SPACECRAFT  MANEUVERS  AND  VIBRATION  ' 


The  momentum  exchange  controller  in  which  the  time  rate  of 
change  of  the  flexible  momentum  relative  to  the  rigid  body  motion 
is  used  as  a  pari  of  the  feedback  control  law  for  maneuvers  and  vi¬ 
bration  suppression  of  flexible  systems  is  introduced.  A  particular 
model  of  a  rigid  hub  (base)  with  a  cantilevered  flexible  appendage 
undergoing  a  single-axis  maneuver  is  considered.  The  feedback 
control  on  the  hub  includes  the  rigid  body  motion  and  the  time  rate 
of  change  of  that  part  of  the  flexible  momentum  resulting  from  flex¬ 
ible  vibration.  The  lower  and  upper  bounds  of  the  control  law  for  sys¬ 
tem  Lyapunov  stability  are  obtained,  and  the  relationship  of  the  con¬ 
trol  law  with  the  energy  Lyapunov  test  function  is  established.  With 
the  presence  of  this  feedback  control,  an  additional  independent 
flexible  control  system  acting  on  the  flexible  parts  can  be  designed 
for  further  vibration  suppression  This  hybrid  control  system  can  be 
applied  to  both  stationkeeping  and  large  angle  maneuvers  Both 
analytical  and  numerical  results  are  presented  to  show  the  theoreti¬ 
cal  and  practical  merit  of  this  approach. 

2.1  Introduction 

The  control  of  flexible  systems  usually  requires  the  controller  to  provide  the  control 
effort  for  maneuvering  or  targeting  of  the  flexible  systems/subsystems  (such  as  anten¬ 
nas  and  telerobols)  with  the  simultaneous  vibration  suppression.  For  instance,  one  of 
the  major  problems  with  robotic  manipulator  systems  is  the  assembly  time  lost  while 
waiting  for  the  vibration  suppression  of  the  end  of  the  manipulator  arm. 


The  mathematical  modelling  of  flexible  systems  itivolves  the  coupling  of  fitiite  and 
infinite  ordered  mathematical  subsystems-the  rigid  body  motions  described  by  ordi¬ 
nary  differential  equations  and  the  flexible  motions  described  by  partial  differential 
equations.The  control  objective  is  to  maneuver  the  rigid  body  and  suppress  the  flexible 
vibration  simultaneously  for  targeting,  or  to  supprass  both  motions  for  stationkeeping 
Most  previous  approaches  to  this  problem  resorted  to  the  discretization  of  the  continu¬ 
ous  systems  into  finite  dimensional  systems  as  the  first  stage  of  the  process  for  control 
law  design  It  is  commonly  recognized  that  the  discretization  procedure  inevitably  in¬ 
volves  modeling  errors  and  control  spillover  issues  Also,  since  the  frequencies  of  the 
flexible  vibration  are  usually  much  higher  than  those  of  rigid  body  motions,  the  coupling 
of  these  two  motions  will  result  in  the  stiffening  of  the  system  differential  equations  As 
a  consequence,  the  numerical  solutions  to  the  traditionally  optimal  approaches,  such 
as  minimizing  a  quadratic  cost  functional  or  minimum  time  bang  bang  control,  are 
usually  difficult  to  solve  because  of  the  inherent  nonlinearity  and  stiffness  of  the  result¬ 
ing  differential  equations.  In  addition,  the  algorithms  and  control  laws  are  usually  sensi¬ 
tive  to  modeling  errors  and  may  lack  robustness. 

In  order  to  suppress  the  flexible  vibrations,  certain  quantities  that  contain  the  informa¬ 
tion  about  the  flexible  motion  of  the  distributed  systems  are  needed  for  the  feedback 
controller.  What  are  these  quantities  for  continuous  mechanical  systems?  From  the 
viewpoint  of  dynamics,  the  rigid  body  motion  and  the  flexible  motion  interact  with  each 
other  by  means  of  momentum  exchange  during  maneuvering;  hence,  that  part  of  the 
flexible  momentum  resulting  from  the  flexible  motion  is  expected  to  contain  the  informa¬ 
tion  about  the  flexible  vibration  and  can  be  used  as  a  part  of  the  feedback  law  The  main 
purpose  of  the  present  paper  is  to  introduce  the  so-called  momentum  exchange  con¬ 
trol  concept  which  is  illustrated  by  the  application  to  a  particular  model  of  a  rigid  hub 
with  a  flexible  appendage  beam  undergoing  a  single-  axis  maneuver.  The  proposed 
control  acting  on  the  hub  includes  the  feedback  of  the  rigid  body  motion  and  the  time 
rate  of  change  of  the  flexible  momenlurn  of  the  beam  (combination  of  the  flexible 
modes)  We  will  strictly  prove  the  stability  of  this  control  law  by  using  a  Lyapunov  test 
function,  and  the  lower  and  upper  bounds  of  the  control  law  for  system  Lyapunov  stabil¬ 
ity  will  be  obtained.  Many  interesting  and  important  implications  of  this  control  concept 
will  be  demonstrated  with  the  development 

Concerning  the  related  research,  attention  should  be  directed  to  recent  impressive 
work  by  Junkins.  et  al'  ,  in  which  the  maneuver  and  vibration  suppression  of  a  rigid 
hub  with  flexible  appendage  beams  was  considered,  and  the  proposed  control  law  in¬ 
volved  the  feedback  of  the  combination  of  the  roof  shear  and  bending  moment  of  the 
beam  which  theoretically  was  equivalent  to  the  time  rate  of  change  of  the  angular  mo¬ 
mentum  of  the  beam  -  rigid  plus  flexible  momentum  The  idea  was  innovative  and  the 
resulting  control  was  easily  implemented  from  the  viewpoint  of  measurement;  however, 
from  the  viewpoint  of  the  present  investigation,  it  appears  that  the  inclusion  of  the  rigid 
momentum  of  the  beam  into  the  feedback  control  law  is  redundant  and  may  corrupt  the 
information  about  the  flexible  vibration  with  possible  degradation  of  the  system  per¬ 
formance,  This  argument  is  further  supported  by  the  fact  that  there  is  no  finite  upper 
bound  limit  ot  the  control  gains  required  for  system  stability  in  Ref  [1]  The  correspond- 


ing  version  of  the  development  of  Junkins'  control  (aw  is  given  in  Appendix  A  for  the 
purpose  of  comparison. 

Another  aim  of  the  present  paper  is  to  desigti  an  additional  independent  flexible  feed¬ 
back  control  system  for  further  vibration  control.  Independent  flexible  feedback  control 
means  that  the  feedback  of  this  control  system  does  not  depend  either  upon  the  open 
loop  rigid  body  maneuver  strategy  or  on  the  precalculated  reference  flexible  motion :  the 
feedback  depends  only  on  the  instantaneous  flexible  displacements  and  velocities  Ba- 
ruh  and  Siverberg  first  suggested  separating  the  maneuver  and  vibration  control  prob¬ 
lem'  .  It  was  proposed  to  suppress  the  vibrations  by  the  natural  modal  control  which  did 
not  affect  the  rigid  body  maneuvers  in  the  sense  that  the  control  forces  conserve  the 
linear  and  angular  momentum  of  the  rigid  body  motion 

Together  with  the  presence  of  the  control  acting  on  the  rigid  hub  in  the  current  paper, 
additional  natural  modal  control  acting  on  the  flexible  beam  for  further  vibration  sup¬ 
pression  will  also  be  employed.  This  hybrid  control  system  is  said  to  be  independent 
in  the  above  sense,  and  it  will  be  shown  that  the  system  stability  is  guaranteed  as  long 
as  the  flexible  control  forces  or  torques  are  energy -dissipating. 

In  summary,  the  objectives  of  this  paper  are  to:  (1)  present  a  kind  of  momentum  ex¬ 
change  feedback  control  concept  and  prove  the  system  stability  by  the  Lyapunov  direct 
method:  (2)  design  an  additional  independent  flexible  control  system  for  further  vibra¬ 
tion  control;  (3)  apply  this  hybrid  control  system  to  the  large  angle  maneuver  problem 
of  the  flexible  system;  and  (4)  demonstrate  the  feasibility  and  practical  merit  of  this  ap¬ 
proach  by  numerical  simulation. 

2 . 2  Momentum  Exchange  Feedback  Control 

A  rigid  hub  (base)  with  a  cantilevered  flexible  appendage  is  considered.  The  append¬ 
age  is  considered  to  be  a  uniform  flexible  beam,  and  the  shear  deformation  and  distrib¬ 
uted  rotary  inertia  are  neglected.  The  beam  is  assumed  to  undergo  small  deformations 
so  that  linear  elastic  theory  can  be  applied.  The  axial  deformation  is  also  neglected.  Mo¬ 
tion  is  restricted  to  the  horizontal  plane,  and.  at  this  stage  of  development,  only  a  control 
torque  acting  on  the  hub  is  considered. 

With  reference  to  Fig.  1 ,  the  position 
vector  of  a  differential  element,  dm.  in 
the  undeformed  body-fixed  frame,  can 
be  expressed  as 

r  =  X  -(-  y  i  (1) 

The  inertial  velocity  of  dm  can  be  obtained 

r  =  -yei^-»-{xe  +  y)[  (2) 
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Therefore,  the  total  inertial  angular  momentum  with  respect  to  point  0  can  be  written 
as; 


H  =  I  ^  9  k  +  J[  X  (  y  +x  0 )  i-y  ■  n  ]  dm  k 

/ 

=  [  (  I,  +  lh  )  e  +  fp(  X  y  +y'  0  )dx  ]  k 


(3) 


where 


I  h  is  the  moment  of  inertia  of  the  rigid  hub  w  r.t.  the  point  O 

/ 

I  b  =  /  n  is  the  moment  of  inertia  of  the  undeforined  beam  . 

/  ^ 

That  part  of  the  angular  momentum  associated  with  the  beam  s  flexible  motion  is  :le- 

I 

noted  H,  =  ;p(xy  +  y’0)dx.  and  is  called  the  flexible  momentum  of  the  beam. 

After  applying  the  Euler-Newtonian  equation  of  motion,  one  can  obtain  the  rotational 
motion  equation  for  the  rigid  body  motion  as; 


( 'h  +'b ) e  1-,^'  =  u(t) 

where  u(t)  is  the  control  torque  acting  on  the  hub. 


(4) 


The  partial  differential  equation  governing  the  elastic  vibration  of  the  beam  can  be  writ¬ 
ten  as; 


P( 


3^ 

3tJ 


-I-  X 


dt*  ^ 


+  HOT 


(5) 


The  H  O  T,  denotes  higher  order  terms  which  include  other  known  effects  [  such  as  rota- 

•  tional  stiffening  (  y  eM  ,  shear  deformation,  etc.).  The  fundamental  development  given 

here  does  not  consider  these  higher  order  terms 


The  boundary  conditions  on  Eq.(5)  are 


• 

at  X  =  /^ 

y(t. 

at  X  =  / 

3^1 

3x*l 

.alyi 

X  -/  3x’  1  ,<=  /  -  0 

^  A  feedback  control  law  for  stationkeeping  is  introduced  as; 

u(t)  = -k,9  -  kjO  -  Y  H,  (7) 

The  purpose  of  the  first  two  feedback  terms  in  (7)  is  straight-forward,  providing  the  feed- 
I  back  control  for  the  rigid  body  motion  of  the  system;  the  last  term  in  (7)  is  the  feedback 

I#  of  the  time  rate  of  change  of  the  flexible  momentum  of  the  beam  and  involves  nonlinear 

output  of  the  flexible  motions. 
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Note  that  the  form  of  this  feedback  control  is  different  from  that  presented  in  Junkins. 
at  el'  ^ ,  where  the  control  involved  the  feedback  of  both  the  root  shear  and  bending  mo¬ 
ment  of  the  beam,  and  the  combination  of  these  moments  was  equal  to  the  time  rate 
of  change  of  the  angular  momentum  of  the  beam  rigid  plus  flexible  momentum 

However,  the  selection  of  the  Lyapunov  function  for  system  stability  here  is  similar 
to  that  in  Ref.[1];  the  slight  difference  lies  in  the  first  term  of  the  Liapunov  function  The 
Liapunov  test  function  ( "error  energy")  is  taken  as  the  following  form; 


The  first  term  is  the  kinetic  energy  of  the  rigid  body  motion  of  the  system:  the  second 

term  “  the  torsional  spring  energy "  is  added  to  make  the  final  state  ( 0  .  ) 

=  (0,0, 0,0)  be  the  global  minimum  of  V;  the  third  term  is  the  kinetic  and  potential  energy 

of  the  beam.  The  weighting  coefficients  a ,  are  included  to  allow  relative  emphasis  upon 
the  three  contributors  to  the  “  error  energy  "  of  the  system 


It  is  obvious  by  inspection  that  imposing  ai>o.  guarantees  that  V  >  o .  and  that  the 
global  minimum  of  V  =  0  occurs  only  at  the  zero  state.  Differentiation,  substitution  of  the 
equations  of  motion  (5)  and  boundary  conditions  (6),  and  some  calculus  manipulations 
lead  to 


V  =  a,(l,  -i-IJ  O9+a,90  +a,9  El  (/„| 
Furthermore,  for  the  rigid  hub.  we  have 


3jy 


a'y 


ay 


Xr-/ 


(9) 


(10) 


Upon  substituting  (4),  (7).  and  (10)  into  (9),  one  can  obtair  s 


v=  {(•a,-k,(  a,-f.-p^  )  I  e  - k^(  +  )g  -(( 3,-7%  )+y(  a,+-^  )]  H ,  }^ 

’  h  '  b 


+  1."''  Ih  +  Ih'"  ‘'^*'"1  -n  '"''"’^L+l,  ''dt 


£2 

h  T'b 


(11) 


In  order  to  meet  the  requirement  that  V  <  0  to  guarantee  stability,  two  of  the  coefficients 
in  (1 1)  are  set  to  zero; 


3,-1- 


U  +  'b 


•  h+l  b 
)+Y(a,+ 


)=  0 


( a,-  \  -Lw/-  .  -£2!^ 


)  =  0 


(12) 

(13) 
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In  this  case,  v  becomes 


aj 


V  =  -kJa +— -j- (ifk.>o) 

I  h  ‘r  I  u 


Eqs. (12-14)  lead  to  the  following  control  gain  requirements; 

k  =  ^  h  +  *b  )  ^ 

3,lh  +a,( l|,  +lf, ) 

_  a, 1 1,  -ajlh  +lb  ) 
a,l^+a,(l|^  +  lf, ) 


(14) 

(15) 

(16) 

(17) 


Since  a,  >  o  ,  it  can  be  easily  proven  that 

<y 

'  b 

and  7  = -1  when  a, =  0;  7=lh/lb  whena,=  o. 


Thus  a  finite  lower  and  upper  bound  re¬ 
stricts  the  size  of  the  control  gain.  v.  It 
should  be  noted  that  there  was  no  such 
finite  upper  bound  limit  in  Ref.  [1.21. 
Fig.2  shows  7  as  a  function  of  a,  /  a, 


Now  attention  is  directed  to  the  significant  physical  meaning  of  this  feedback  control 
law.  According  to  Eq.(8).  a,  /a,  corresponds  to  the  relative  contribution  of  the  flexible 
motion  and  the  rigid  body  motion  to  the  system  error  energy.  When  a ,  =0  which  corre¬ 
sponds  to  7  =-l  .  the  “  error  energy  ”  Lyapunov  function  does  not  contain  the  flexible 
energy  term.  In  other  words,  the  emphasis  is  totally  put  on  the  control  of  the  rigid  body 
motion.  For  the  case  where  v=~  1  it  can  be  seen  from  Eqs  (4)  and  (7)  that  the  effect  of 
the  flexible  momentum,  H, ,  is  removed  from  the  closed -loop  rigid  body  motion  equa¬ 
tion.  Thus,  the  controller  completely  compensates  for  the  flexible  motion  and  the  rigid 
body  motion  is  not  disturbed  by  the  flexible  vibration.  Here  a  two-lumped  mass  beam 
model  is  used  to  approximate  the  continuous  beam  (Fig  3);  however,  the  theoretical  de¬ 
velopment  and  the  validation  of  the  results  do  not  depend  on  the  discretization  proce¬ 
dure  Fig. 4  and  Fig.  5  show  the  transient  responses  of  the  rigid  body  motion  and  the  flex¬ 
ible  motion,  respectively,  wherein  the  point  2  is  the  outermost  lumped  mass.  The  three 
cases  considered  are:  v=- 1 ,  maximum  emphasis  on  the  rigid  body  motion  (complete 
compensation  for  the  flexible  motion);  v=n.  no  feedback  of  the  flexible  momentum; 
i’=lh/lb=I  .83 .  maximum  emphasis  on  the  flexible  motion.  It  can  be  seen  that  with  the 
increase  of  v  (larger  relative  emphasis  on  the  suppression  of  the  flexible  motion)  the 
dissipation  of  the  flexible  vibrations  is  better,  and  that  the  rigid  body  motion  is  subjected 


(18) 


V 


Fig.  2  V  as  a  function  of  a, /a, 


to  more  disturbance,  but  the  disturbance  is  not  so  large  as  to  degrade  the  rigid  body 
motion  seriously.  Therefore,  in  practical  applications  the  control  gain  can  be  selected 
according  to  the  requirement  for  the  rigid  body  motion  and  the  flexible  vibration  sup 
pression.  and  it  is  a  kind  of  trade-  off.  The  selection  of  the  control  gain  can  be  further 
enhanced  by  combining  some  optimum  <;rlteria‘ 

Also,  it  should  be  noted  that  bounded  restrictions  for  the  gain  r  are  only  sufficient 
conditions  for  system  stability;  however,  from  the  present  simulation,  it  is  seen  that  these 
restrictions  are  close  to  the  necessary  conditions  for  system  stability,  especially  the  low¬ 
er  bound  condition.  Fig.  6  and  Fig.  7  show  that  when  the  v  exceeds  the  lower  bound  the 
flexible  motion  becomes  unstable  and  vice  versa,  when  the  v  exceeds  the  upper  bound 
the  rigid  body  motion  diverges 

The  system  performance  will  generally  depend  on  the  ratio  1 1,  / 1  b .  Fig.  8  and  Fig 
9  show  the  system  responses  corresponding  to  three  different  ratios.  It  is  seen  that  if 
this  ratio  increases,  the  flexible  vibration  will  be  usually  easier  to  be  suppressed  and  the 
rigid  body  motion  will  also  decay  more  quickly  with  a  little  bit  larger  overshoot  at  the  be¬ 
ginning.  Also,  from  the  simulation  results,  it  is  seen  that  this  control  law  is  robust  in  the 
presence  of  other  unmodelled  factors  and  modelling  errors,  and  the  robustness  will  de¬ 
grade  when  the  control  gain  i»  approaches  its  boundaries. 

2 . 3  independent  Flexible  Control  for  Elastic  Vibration 

In  Section  II,  it  has  been  demonstrated  that  the  feedback  controller  located  on  the 
rigid  part  of  the  spacecraft  can  perform  the  maneuver  and  simultaneous  vibration  sup¬ 
pression,  and  the  system  stability  is  guaranteed  as  long  as  the  control  law  gains  are 
within  a  certain  range.  However,  it  is  of  interest  to  further  control  and  suppress  the  flex¬ 
ible  vibration  while  the  spacecraft  is  going  through  large  angle  rotations  and  large  tr  ans- 
lations.  To  achieve  this,  another  independent  flexible  control  system  acting  on  the  flex¬ 
ible  appendage  is  proposed. 

Baruh  and  Silverberg^  proposed  a  natural  modal  control  for  the  independent  elastic 
control  system.  They  have  proven  that  the  natural  control  forces  conserve  the  linear  and 
angular  momentum  of  the  spacecraft  and.  therefore,  do  not  appear  on  the  right-hand 
side  of  the  rigid  body  equations.  However,  the  rigid  body  motion  is  still  affected  by  the 
flexible  motion  via  the  kinematic  and  dynamic  coupling  which  is  sometimes  very  impor¬ 
tant.  With  the  feedback  control  presented  in  Section  II,  it  will  be  proved  that  the  system 
stability  is  guaranteed  as  long  as  the  elastic  modal  control  system  is  an  energy -dissi¬ 
pating  system.  To  prove  this,  it  is  assumed  that  u^  is  the  feedback  natural  modal  control 
on  the  flexible  appendage;  hence,  can  be  written  as 

Ue=Irt>rfr  (19) 

where  (t>r  is  the  rth  natural  modal  shape  corresponding  to  the  flexible  motion  and  f , 
is  the  corresponding  modal  control  force  which  can  be  designed  in  the  modal  space 
using  various  methods,  such  as  optimal  LQR,  pole  assignment,  or  velocity  control 
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(20) 


In  this  case.  Eq.(5)  becomes 

p  (  —  +  X  —  >  El  -  U 
'  '  3t=  dt-  dX' 

The  Lyapunov  function  is  taken  the  same  as  (9)  Using  the  orthogonality  relation  of  the 
natural  modal  control  (see  details  in  Ref  [3]).  and  after  similar  manipulations,  the  deriva¬ 
tive  of  V  can  be  obtained  as 

/ 

V  = '■'Qht-hand  side  of  Eq.(14)  -i-  J'pyu^dx  (21) 


Therefore,  if 


f 

J py  u„dx  <  0 


which  means  the  control  forces  are  energy-dissipating,  then  V  <  0  and  ihc  ^vrtcin  i<; 
stable. 

Note  that  truncation  of  natural  modes  (finite  control)  does  not  change  the  discussion 
here.  However,  in  reality  it  is  almost  impossible  to  obtain  the  exact  natural  mode  shapes 
which  are  usually  approximated  by  discretizing  the  continuous  system  model;  the  pro¬ 
cess  of  discretization  inevitably  introduces  natural  control  errors.  From  the  simulation 
results  here,  the  momentum  exchange  control  law  is  robust  and  can  endure  these  er¬ 
rors  induced  from  approximating  the  natural  modal  shapes.  The  simulation  results 
show  that  ttie  addition  of  this  elastic  control  system  further  suppresses  the  elastic  vibra¬ 
tion  and  improves  the  system  responses  (Figs.  10  -1 1). 

2 . 4  Large  Angle  Maneuvers 

For  large  angle  maneuvers,  a  reference  open-loop  rigid  body  maneuver  is  adopted, 
denoted  by  the  subscripts  as  {QrC’  .  ;  ^  satisfying 

(  1^,-t-lb  )  9,  =  uut)  (23) 

and  the  prescribed  boundary  conditions.  The  variables  without  subscripts  represent 
the  actual  solutions  to  the  system  equations  (4),  (7).  and  (20). 

The  tracking  control  law  takes  the  following  form; 

U(t)  =  Ur  -k,(e-0r)-k^(e-er)  - 

The  "error  energy”  Lyapunov  function  is  modified  as  ^ 


■fB&dx  1  } 
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The  same  manipulations  as  in  Sect.  II  lead  to 


v  =  {[  a,-k,(  a ^  )  ]  (ft-e,)  -  l<,(  a,+  )(n-o.) 

L  +  I  h  1 1.  'f  '  1 1, 


3  I  3  I 

-((a-r^-  )+'^'(  }{A-A,)+a,  f ny  u^dx  -aj pxydx  a. 

^  h  h  h 


Again,  the  tracking  control  laws  are  subjected  to  the  same  requirements  as  in  ( 1 5).  ( 1 6). 
and  (17).  In  this  case, 

a  I  •  2  r  r 

V  =-k,(a,+  -j-^  )(e-A.)  +a,  Jpy  u,dx  a,Jpxydx  a,  ^27) 

in  which  the  first  term  is  semi-negative  definite:  the  second  term  is  negative  definite  if 
the  flexible  control  force  is  energy-dissipating;  and  the  sign  of  the  third  term  is  generally 

varying.  Therefore,  as  a  whole,  we  cannot  say  the  \  is  negative.  In  effect,  the  system 
has  continuous  inputs  and  the  flexible  vibration  will  be  excited  during  large  angle  ma¬ 
neuvers.  However,  it  is  hoped  that  with  this  hybrid  control  the  flexible  vibration  will  be 
suppressed  and  that  bounded-input/bounded-output  viewpoint  of  stability  can  be  es¬ 
tablished. 

It  is  not  the  purpose  to  discuss  the  open -loop  maneuver  strategy  here.  A  reveiew 
about  this  subject  can  be  found  in  Ref.(4).  For  the  purpose  of  demonstration,  a  mini- 
mum-time  nianeuver  strategy  which  is  a  bang -bang  law  is  considered.  It  is  recognized 
that  a  minimum-time  maneuver  bang-  bang  law  will  excite  the  flexible  vibrations  exces¬ 
sively.  and  it  is  intentionally  used  here  to  show  the  effectiveness  of  the  proposed  hybrid 
feedback  control 

For  quick  large  angle  maneuvers,  k,  .  k^  are  usually  taken  to  be  zero.  Notice  that  in 
this  case  the  remaining  feedback  term  ( -I’^j  H , )  does  not  depend  on  the  open-loop 
maneuver  strategy  and  the  precalculated  corresponding  reference  flexible  motion.  The 
results  of  a  90  degree  minimum-time  maneuver  are  shown  in  Figs.  1 2- 1 4.  It  is  apparent 
from  Fig.  14  that  the  deflection  of  the  outermost  lumped  mass  (point  2)  is  reduced  by 
almost  half  when  the  flexible  momentum  is  included  as  a  feedback  term  in  the  tracking 
control  law. 


2 . 5  Conclusion 

A  kind  of  momentum  exchange  control  law  for  maneuvers  and  vibration  suppression 
of  flexible  systems,  wherein  the  time  rate  of  change  of  the  flexible  momentum  is  used 
as  a  part  of  the  feedback  control  law,  has  been  presented.  The  theoretical  analysis  for 
system  stability  by  using  a  Lyapunov  function  is  rigorous  without  resort  to  the  disctreti- 
zation  of  the  continuous  system,  and  the  lower  and  upper  bounds  of  the  control  gains 
for  system  stability  have  been  obtained.  For  further  vibration  suppression,  the  natural 


modal  flexible  control  system  is  used  in  conjunction  with  the  momentum  exchange  con¬ 
troller.  This  hybrid  control  system  has  also  been  applied  to  the  large  angle  maneuver 
of  the  flexible  system  in  which  the  minimum  time  open  loop  maneuver  strategy  is 
used.  Simulation  results  have  shown  the  effectiveness  of  this  control  system  for  both 
stationkeeping  and  minimum-time  maneuvers 

This  study  is  interesting  and  significant  for  understanding  the  behavior  and  control 
of  distributed  mechanical  systems.  It  is  noteworthy  to  see  that  just  the  single  output 
feedback  of  the  time  rate  of  change  of  the  flexible  momentum  contains  enough  informa¬ 
tion  for  the  controller  to  suppress  the  vibration  of  continuous  systems  (infinite  order  sys¬ 
tems).  Also,  the  momentum  exchange  control  law  tends  to  be  of  a  nonlinear  form  of  the 
flexible  displacements.  This  study  has  demonstrated  the  prospect  and  efficacy  of  using 
nonlinear  control  in  some  specific  systems  such  as  mechanical  systems.  Another  key 
point  is  the  usefulness  of  Lyapunov  functions  in  the  control  problems  of  distributed  sys¬ 
tems. 
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Apendix  A.  A  Comparison  of  Junkins'  Control  Law  with  the  Present  Law 


In  Ref.  { 1 ,2],  Junkins.  at  el.  presented  a  Liapunov  control  law  design  for  the  maneuver 
and  vibration  suppression  of  flexible  systems  The  model  considered  was  a  rigid  hub 
with  flexible  appendage  beams  undergoing  a  single  .axis  maneuver,  and  each  beam 
had  a  finite  tip  mass.  For  easy  comparison,  only  one  appendage  and  no  tip  masses  are 
adopted  in  the  following  development  of  Junkins  control  law.  In  terms  of  the  notations 
of  Ref  [1].  the  system  equations  of  motion  ate: 

*h  ^  (28) 

-(  Mo -So/o)  I7,  +  X  ^  (29) 

/, 

+  '  0  "HOT  (30) 

where  (  Mo  So)  denotes  the  bending  moment  and  the  shear  force,  at  the  root  of  the 
beam,  respectively,  and  (30)  satisfied  the  same  boundary  conditions  as  (6) 

After  substitution  of  (29)  into  (28).  one  can  obtain 

dh  )  9  +  =  U(t)  (31) 

H ,  is  the  flexible  momentum  of  the  beam 

The  control  law  of  Ref.(1]  is  taken  as  the  form: 

u= -[  g,  9 +9^9  Mo -S./J  1  (32) 

After  substitution  of  (29)  into  (32),  the  control  law  can  be  rewritten  as 

U  =  -(  g,  9+g29  +  g3(  l^9  +  |-^H)  )]  (33) 

Notice  this  control  law  is  different  from  (7)  in  that  the  third  term  of  (33)  involves  the  extra 
time  rate  of  change  of  the  rigid  momentum  of  the  beam 

In  order  to  show  the  difference,  we  use  the  same  Lyapunov  function  (8)  as  the  test 
function.  In  this  case,  the  same  manipulations  as  in  Sect  II  lead  to  the  following  require¬ 
ments  for  the  control  gains: 


g 


92  >0 

_  t^g3)M 

3,1  +a,  ( I H  + 1 ) 


(34) 


In  this  case,  the  derivative  of  the  Lyapunov  function 


V 


a.dn  -Hb)  2 


(35) 


From  (34), 


93  >  - 1  and  without  finite  upper  bound 


Also,  from  (34)  and  (35).  it  is  seen  that  if  a,  is  fixed  and  if  a,  is  increased  (  means  the 
increase  of  the  relative  contribution  of  the  flexible  motion  to  the  error  energy  Lyapunov 
function),  the  gain  g,  will  become  large;  however,  the  v  is  unchanged  This  means 
physically  that  increasing  the  control  gain  g, ,  unlike  increasing  the  control  gain  v.  will 
not  automatically  increase  the  suppression  of  the  flexible  motion.  Therefore,  it  is  con¬ 
cluded  that  the  inclusion  of  the  rigid  momentum  into  the  feedback  control  law  is  not  nec¬ 
essary  and  may  corrupt  the  feedback  information  with  possible  degradation  of  the  sys¬ 
tem  performance. 

It  Is  easily  seen  from  Eg.  (29)  that  the  combination  of  the  root  shear  and  bending 
moment  is  equal  to  the  time  rate  of  change  of  the  rigid  momentum  plus  the  flexible  mo¬ 
mentum  of  the  beam.  Therefore,  the  realization  of  the  current  proposed  control  law  is 
not  difficult;  the  flexible  momentum  feedback  term  can  be  realized  by  measuring  the 
root  shear  force,  root  bending  moment,  and  the  angular  acceleration  of  the  hub. 
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Table  1  Spacecraft  Dimensions,  Appendage  Material 


Radius  of  the  rigid  central  body,  / « 

1  00  m 

Length  of  the  appendage.  / 

20.00  m 

Appendage  material  stiffness.  El 

1500.00  N-m  ^ 

Appendage  material  density,  p 

0.04096  rn 

Mass  of  the  rigid  hub 

400.00  kg 

Total  rotational  inertia,  Ih+lb 

309.1863  kg-m^ 

Inertia  ratio.  Ih/lj, 

1.83 

Fig.  3 


Lumped-Mass  Model  for  Numerical  Simulations 


2.13 


Fig.  4  Rigid  Body  Rotation  Responses  in  Stationkeeping 
under  Different  Stable  Control  Laws 


Fig.  5  Flexible  Vibration  Responses  in  Stationkeeping  under 
Different  Stable  Control  Laws  (IC’s  the  same  as  in  Fig. 4) 


Note;  k,  =  618and  k2  =  773  in  Figs.  4-11  and  k,  .  k-  =0  in  Figs  12  14 
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Deflection  of  Point  2  (ml 


Fig.  6  Flexible  Vibration  Response  when  i-  Exceeds  the  Lower 


Bound  of  Stability  Range  (IC’s  the  same  as  in  Fig  4) 


•  Upper  Boundof  Stability  Range  (IC  ’s  the  same  as  in  Fig  4) 
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Deflection  of  Point  2  (m) 


Time  is) 

Fig.  8  Transient  Responses  when  the  Ratio  Changes 


(IC’s  the  same  as  in  Fig.^) 
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_  Deflection  of  Point  2  (m) 

Ll,(0)  =  2  m 

2  -  e  (0)  =  o’ 


i’=1,83 

without  the  natural  modal  control 
with  the  natural  modal  control 


2  4  6  8  10  12  14  16  18  20 

Time  (8) 

Fig.  10  Transient  Responses  with  and  without  the  Natural 

Modal  Control 


Theta  (deg.) 


i'=1.83 

. without  the  natural  modal  control 

—  with  the  natural  modal  control 


2  4  6  81012141618 

Time  (a) 

Fig.  1 1  Transient  Responses  with  and  without  the  Natural 
Modal  Control  (IC's  the  sme  as  in  Fig.  10) 
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3.  Adaptive  Feedback  Control  for  the  Retrieval  of  an  Orbiting 
Tethered  Antenna/Reflector  System 


ABSTRACT 

An  adaptive  feedback  control  concept  for  the  retrieval  of  an  orbiting  tether 
antenna/reflector  system  is  introduced.  In  conjunction  with  an  open-loop 
exponentially  decreasing  tether  commanded  length  control,  an  attitude  control 
system  acting  on  the  antenna  is  used  for  stabilizing  the  antenna  as  well  as 
suppressing  the  tether  swing  librationsthrough  the  control  of  the  attitude  motion  of  the 
antenna(base).  The  momentum  exchange  control  technique  wherein  the  time  rate  of 
change  of  the  part  of  the  angular  momentum  associated  with  the  relative  motion  of  the 
tether  is  used  as  a  part  of  the  feedback  is  applied  to  the  control  law  design.  The 
selection  of  the  dimensionless  control  gains  is  based  both  on  the  consideration  of  the 
stability  of  the  system  linearized  about  a  reference  length,  and  on  the  numerical 
simulation  of  the  actual  nonlinear,  time-varying  system.  The  control  gains  are  adapted 
to  the  instantaneously  changing  tether  length.  Numerical  results  are  used  to  verify  the 
control  concept. 


3.1  INTRODUCTION 

Since  the  early  1970's  a  number  of  very  large  space  antennas  have  been 
proposed  for  power  transmission,  astronautical  research,  and  communications.  The 
gravity  stabilized  configuration  is  particularly  suited  for  very  large  flexible  systems  to 
alleviate  the  problems  associated  with  the  active  attitude  control  of  very  large 
structures.  An  orbiting  tethered  antenna/reflector  system  has  been  proposed  for 
possible  communication  or  military  use[1,2].  The  subsatellite  is  nominally  deployed 
below  the  antenna  along  the  yaw  axis  at  a  sufficient  distance  (1km)  to  provide  a 
favorable  composite  moment  of  inertia  ratio  for  gravitational  stabilization.  It  is 
assumed  that  the  tether  would  be  deployed  through  the  end  of  a  rigid  boom  attached 
to  the  shell(antenna)  apex.  The  purpose  of  the  present  paper  is  to  investigate  the 
control  strategy  for  the  retrieval  of  the  tethered  antenna  system. 

It  is  well-known  that  retrieval  is  basically  unstable  because  of  the  excitation  of 
both  in-plane  and  out-of-plane  librations  and  insufficient  tension  during  the  terminal 


phase.  Many  investigators  have  studied  the  tether  deployment/retrieval  problem 
Most  of  the  previous  approaches  on  letrieval  were  based  on  modulation  of  the  tether 
tension:  however,  tension  control  will  not  be  effective  for  short  tethers  because  of  the 
low  level  of  the  gravity  gradient,  and  it  is  found  that  the  tension  control  without  a 
constraint  on  the  minimum  tension  level  could  easily  lead  to  a  slack  tether. 

Based  on  the  momentum  exchange  control  concept  introduced  in  a  recent 
paper[3],  a  new  control  strategy  for  suppression  of  the  tether  swing  librations  during 
retrieval  is  proposed  in  the  present  paper.  It  is  proposed  that  together  with  an 
open-loop  exponentially  decreasing  tether  commanded  length  control  an  attitude 
feedback  control  acting  on  the  antenna  be  utilized  for  stabilizing  the  antenna  as  well 
as  suppressing  the  tether  in-plane  swing  libration.  The  resulting  feedback  control  law 
would  include  two  components.  One  component  involves  the  linear  feedback  of  the 
antenna  rigid  body  motions  for  the  purpose  of  guaranteeing  the  stability  and 
convergence  of  the  antenna  attitude  motions.  The  dimensionless  feedback  gains  of 
this  part  are  adapted  to  the  initial  and  final  phases  of  the  retrieval.  The  other 
component  involves  a  feedback  of  the  time  rate  of  change  of  the  part  of  the  angular 
momentum  associated  with  the  in-plane  motion  of  the  tether  relative  to  the  local 
vertical.  The  purpose  of  this  part  of  the  feedback  law  is  to  provide  the  information 
about  the  tether  in-plane  swing  motion  for  the  controller  to  suppress  the  tether  swing 
motion  (momentum  exchange  control).  The  dimensionless  gain  for  this  part  is 
selected  based  both  on  the  consideration  of  the  stability  of  the  system  linearized  about 
a  reference  length  and  on  the  numerical  simulation  of  the  actual  nonlinear, 
time-varying  system.  The  gain  is  adapted  to  the  instantaneously  changing  tether 
reference  length.  The  possible  practical  advantage  of  this  new  approach  over 
traditional  tension  control  during  retrieval  is  that  the  longitudinal  vibration  induced  by 
this  control  method  should  be  less  than  that  induced  by  tension  control;  in  general, 
tension  control  will  directly  excite  the  longitudinal  vibration.  It  has  been  observed  that 
longitudinal  vibrations  and  wave  propagation  may  present  serious  problems  for  tether 
operations,  and  the  solution  to  these  problems  is  being  sought. 


THE  SYSTEM  EQUATIONS  OF  MOTION 

For  system  modelling  the  following  assumptions  are  made; 

•  The  shift  of  the  center  of  mass  of  the  system  is  neglected.  The  center  of  mass  is 
assumed  to  lie  at  the  same  point  as  the  center  of  mass  of  the  antenna,  and  the  center 
of  mass  of  the  system  is  moving  in  a  circular  orbit. 

•  The  flexibility  of  the  shell,  boom,  and  tether  is  neglected  and  the  tether  is  assumed 
uniform  with  a  constant  mass  per  length. 

•  The  symmetry  axis  of  the  shell  is  nominally  along  the  local  vertical. 

•  The  subsatellite  is  considered  as  a  point  mass. 

The  Newton-Euler  method  has  been  adopted  to  develop  the  system  equations 
of  motion.  The  coordinate  systems  used  in  the  development  of  the  system  equations 


of  motion  are  shown  in  Fig. 2 .  OpX^YoZ,  is  an  orbit-fixed  reference  frame  centered  at 
the  center  of  mass  of  the  shell,  Op,  with  Op  along  the  local  vertical  and  Op  Y,,  along 
the  orbit  normal  opposite  to  the  orbital  angular  velocity  vector.  OpXpYpZp  is  a 
shell-fixed  reference  frame,  Rp  ,  where  OpXp ,  OpYp ,  Op  Zp  are  the  principal  axes  of 
the  shell.  OXYZ  is  the  tether-fixed  reference  frame,  Rt  ,  with  OX  along  the  tether  line, 
where  O  is  the  point  from  which  the  tether  is  deploying.  The  coordinates  of  point  O  in 
the  shell  frame,  Rp ,  are  (  h^ ,  0, 0 ) 

The  Euler  angles,  \1;,  0,  <t),  are  the  yaw,  pitch  ,  and  roll  angles  of  the  shell, 
respectively,  a.  y  are  in-plane  and  out-of-plane  swing  angles  of  the  tether. 
Therefore,  the  whole  system  has  the  following  degrees  of  freedom; 


\i/,  0,  cf) - rigid  body  motion  of  the  shell 

rt.  y  - librational  motion  of  the  subsatellite 

I  - length  of  the  tether 


The  transformations  from  OpXpYpZp  to  OpXpYpZp  and  from  OpXpYpZp  to  OXYZ  are 
given  by 
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Consider  an  elemental  mass,  dm.  whose  instantaneous  position  vector  from  the 
center  of  the  shell.  Op ,  is  T  (Fig.  2).  The  equation  of  motion  for  dm  can  be  written  as 

a  dm  =  f  dm  +  e  dm  (3) 


3.3 


where 


a  =  inertial  accelei  ation  of  dm 
T  =  gravitational  force  per  unit  mass 
e  =  external  forces  acting  per  unit  mass 
The  gravity  force  in  the  shell  frame,  Rp,  is  given  by  [4] 

r=r  +  MF  (4) 


where  is  the  gravity  force  at  Op  expressed  in  the  frame,  Rp.  and 
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(5) 


where  is  the  orbital  angular  velocity. 

The  vector  equation,  (3) ,  can  be  written  in  the  frame,  Rr ,  as 

f  ^r.“r„+r+2  cox7  X ( oixf )  +  w X  r  -M  r  ]  dm -c dm  =  0  (6) 

where 
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The  shift  of  the  center  of  the  mass  of  the  entire  system  from  the  center  of  the  antenna  is 
neglected,  that  is,  .  After  projecting  equation  (6)  on  the  tether  frame,  R, ,  the 

following  is  obtained 

T[  r  +2  wx?  +w  x(  wxT)  +  w  x  T-M?  ]dm-c|  dm  =  0  (8) 

where  |i  indicates  the  projection  onto  the  frame,  R,  . 

The  expression  for  the  F  of  the  tether  system  is  different  from  that  of  the  shell  due  to 
the  relative  motion  of  the  tether,  so  we  consider  the  tethered  subsatellite  system  and 
the  shell  system,  separately. 

TETHERED  SUBSATELLITE  SYSTEM 

r  =  T'f,  +  h  (9) 
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where  =(  x,  o,  o  )  is  the  position  vector  of  dm  from  O  projected  onto  the  frame. 
R,.  H=(hx,  0.0)  is  the  position  vector  of  point  O  from  the  shell's  center  of  mass.  Or 

Hence, 


7  =  r‘r,+  r'r, 
r=T  'r,+  2t  'r,  +  t  'r, 


(10) 


According  to  vector  algebra 
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After  substitution  of  equations  (10)-(1 1)  into  equation  (8),  there  results 

a,drn--c|  drn  =  0  (12) 

where 
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T[0|h  =  h.. 
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The  expressions  for  some  terms  are  listed  in  Appendix  A. 

Also,  let  I^=J'xdm  H„=J'x-(Jm 

By  projecting  Eq.(12)  along  the  X,  Y,  and  Z  axes  and  integrating  the  result  over  the 
subsatellite  and  tether,  the  translational  equations  of  motion  can  be  obtained  as 

(m,  +  m,)/  +  [,co,,+ (m,+ni,)h,;gi=  F„ 

2  (m,+m,)/  Q^+ +  g,  =  F^.  (14) 

-2(m,  +  m,)/  n,+  I,to„  +  (m,  +  m,)h,  g,= 


where  F,^ ,  F,,, ,  F,^  are  three  components  of  the  force  exerted  on  the  tether  by  the 
shell.  No  other  external  forces  acting  on  the  tether  and  subsatellite  are  assumed.  The 
equations  of  rotational  motion  of  the  tether  can  be  obtained  by  the  following  operation; 


/  rx[Eq.(I2)|  =  0  (15) 


By  projecting  Eq.(1 5)  along  the  Y  and  Z  axes,  respectively,  the  rotational  equations  for 
"the  in-plane  and  out-of-plane  motions  are  obtained  as 


2l,nv/-H„w,, -h,l,g3=  0  (16) 

2Un,/+H^w,, +hJxg>  =  0  (17) 

SHELL  SYSTEM 

The  torques  exerted  on  the  shell  are 

L=L,, +L,  (18) 

where  Lr,  is  the  gravitational  torque  and  L,  the  torque  by  the  tether  whose 
components  in  the  shell  body  frame,  Rp  ,  can  be  obtained  as 
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(20) 


where  J,.vv  are  the  shell’s  principal  moments  of  inertia.  Upon  applying  the 
Euler-Newtonian  equations,  the  equations  of  motion  of  the  shell  attitude  motion  can 
be  obtained  as 

-(  -J,  )Wv  Wj  =  3w“  .s9c0S(t)(  -J, )  +  Lex 

J,,  u),,  _(  _jJa),  oj^  =  -3co;  .secec4)(  J,  -JJ  +  h,(.srtcy  F„.  -srt.syFiv  -ca  F„  )  +  Lev 

h.  CO7.  -( J, -J,  )aJx(0v=  3w;  c"e‘^4)C(t)(  J,  -J,. )  -hx(  '^YF.x  +  c7  Fn. )+  Ler  (21) 


where  the  forces  Fix.F,v.F''z  already  been  obtained  in  Eq.(14)  and 

L  e't  ,  L  ey  ,  L  e?  are  other  external  torques  including  the  control  torques.  The  system 
nonlinear  equations  of  motion  (16),  (17),  and  (21)  are  used  to  obtain  the  system 
linearized  equations  and  to  simulate  the  closed-loop  system. 

SYSTEM  QUASI -LINEARIZED  EQUATIONS  OF  MOTION 


In  the  case  of  retrieval  the  time  rate  of  change  of  the  tether  length  and  its  rate  can 
be  fast  and.  therefore,  are  treated  as  zero  order  terms  (playing  the  same  role  as  the 
constant  terms  in  the  stationkeeping  case)  in  obtaining  the  system  quasi-linearized 
equations  of  motion.  The  resulting  quasi-linearized  nondimensional  system 
equations  can  be  obtained  as 


in-plane 


r  k,  a’*+9  ”+3a  +39  +k  ,(«’+  9'-  i)y-  =  () 

^  ^Q,yQ  +  2k^(a'+Q'-\)Y  +k5«y  =  Uq 


(22) 


k,7”-4)”-44)  +(3+k,)7+k,(7’-(t)'-  vj/)  j-’  =  0 

out-of-plane vj/”-nxvi/-(l+nx)<t)’  =  Li,;, 

k,7"-4)”+k,7-4n,4>-(  l-n,)\l/’+2k^(7  (t>’-'l')-|-  +  kr,Yy-  = 

where  ()’  =  -^()  .  T  =  Wct  ■  ^0  AF  .  and  Lf,f,  are  pitch,  yaw,  and  roll  axis  control 
torques,  respectively,  and  no  other  external  disturbance  torques  are  assumed.  The 
nondimensional  parameters  in  Eqs.  (22)  and  (23)  are: 
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k,  =- 


H, 


H„+hL 


k,= 


•llx 


J* 


IKx+hIv 


(  h=  -  h,  X)  ) 


j*-j; 


j;-j- 


I,  _  hni,,/  .  .  _  hnki  .  rx  _  „  _  .-^x 

k.-— .  K^-  ,  ii,  ,  u, 

Jy  v»/ 

j;  =  j,  ;  j;  =Jv+ hl,+m,,  h-  ;  j;=J,+ hl,  +  m.,  ir 


n 

j; 


It  can  be  seen  from  the  quasi-linearized  equations  that  in  the  case  of  retrieval,  j-  <o  ■ 
the  in-plane  and  out-of-plane  motions  become  unstable  due  to  the  divergent  terms, 
such  as  (x'j-  ,  Q'j-  ,etc.. 


3 . 3  MOMENTUM  EXCHANGE  FEEDBACK  CONTROL  FOR  SYSTEM  IN-PLANE 
MOTION  AT  A  FIXED  TETHER  REFERENCE  LENGTH 

In  Ref.[3],  Li  and  Bainum  have  introduced  a  new  control  concept  called 
momentum  exchange  feedback  control.  The  same  control  concept  is  applied  to  the 
current  problem;  that  is,  the  time  rate  of  change  of  that  part  of  the  angular  momentum 
due  to  the  tether  deviation  motion  from  the  local  vertical  is  used  as  a  part  of  the 
feedback  control  law  for  suppressing  the  tether  in-plane  swing  motion. 

Let  us  first  obtain  the  expression  cf  the  time  rate  of  change  of  the  tether  angular 
momentum.  From  Eq.(20),  the  pitch  torque  exerted  on  the  tether  by  the  shell  is 

M,=  -l\(  saoF  V  -  Srt.sYF,,  -cu  F„  )  (24) 

After  Eq.(24)  is  linearized  about  a  fixed  tether  reference  length,  we  obtain 
M,  =  -h,  L  ( «■>  9”)  t-nv  h;  e”-.^lix(  K )  9  (25) 

Since  (a  +  9)  is  the  angle  of  the  tether  line  from  the  local  vertical  in  the  linear  sense, 
according  to  Newton’s  law,  the  linearized  time  rate  of  change  of  that  part  of  the  angular 
momentum  due  to  the  tether  in-plane  motion  relative  to  the  local  vertical  is 

-K\Aa"+  e”) 

Therefore,  applying  the  momentum  exchange  control  concept,  the  nondimensional 
feedback  control  for  the  in-plane  motion  is  introduced  as 

Ufl  =  -g,  9  -  g,  9  -  V(  k  a"  +  k  ,  9"  )  (26) 

The  first  two  terms  in  Eq.(26)  are  the  linear  feedback  of  the  shell  pitch  angle  and  its 
rate;  the  last  term  is  the  feedback  of  the  linearized  time  rate  of  change  of  the  angular 
momentum  due  to  the  tether  in-plane  motion  relative  to  the  local  vertical,  for  the 
purpose  of  exchanging  the  tether  momentum  with  the  shell  and  absorbing  the  tether 
swing  motion  energy  by  the  control  system. 

The  closed-loop  system  equations  of  the  in-plane  motion  about  a  fixed  tether 
reference  length  become 
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k,  a"+ e"+3  (v  +  3P>  =0 


k,  a"+  e”-  ^  n,  ■-)  =  -  g, 0  -  g .  P‘  -V(  k  , a”  +  k  ,  e"  ) 

Note  that  the  control  goal  is  to  suppress  the  angle  0  and  (ex  +  9).  After  substitution  of 
a  =  3-9  into  Eq.(27),  one  obtains 

k,3’ +(  l-k,)9  ■•■33  =  9 

(l  +  V)k,3”+f  g,-3n,)9=  0 

The  characteristic  equation  of  (28)  is 

f(\)  =  [k,-k,-k,,(  l-k,)T>l  +k,g,v’+|  3(  l-k,)+ ( g,-3nv )k ,  |\- 

+3g,\+3(gi-3n,)  (29) 

The  Routh-Hurwitz  criteria  is  applied  and  leads  to  following  requirements  for  the 
control  gains: 

g,  >  30,  (30) 

(31) 

-  1  <  V  <  -  (32) 


kdl-k,) 

3.4  RETRIEVAL 

IN-PLANE  MOTION  CONTROL 

The  linear  feedback  control  law  (26)  is  applied  to  the  retrieval  where  the  system 
involves  complex  nonlinear  and  time-variant  ordinary  differential  equations.  However, 
it  is  expected  that  for  a  fairly  moderate  retrieval  rate,  the  closed-loop  system 
responses  will  be  bounded  when  the  control  gain,  v.  is  within  some  range  as  indicated 
by  (32).  Numerical  simulations  are  used  to  select  the  control  gain  v  which  will  generate 
bounded  or  convergent  system  responses.  The  system  simulation  parameters  are 
shown  in  Table  1 . 

TABLE  1.  SYSTEM  SIMULATION  PARAMETERS 


m, 

loooo  kg 

500  kg 

8.35  kg 

h 

a 

Wc 

0.08  km 

0.1  km 

0.001  178  rad/s 

The  decreasing  commanded  length  is  taken  as  the  exponential  form: 


/  =  /oC 


In  =  I  km 
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In  general,  the  range  of  the  control  gain  r  for  bounded  system  responses  will 
depend  upon  the  retrieval  rate  c  .  For  a  particular  retrieval  rate  c  =  3  45xl5’(s  ' )  for 
example,  it  is  observed  (based  on  numerical  simulations)  that  the  system  is  bounded 
when  V'  lies  within  the  interval 


1  <  v< 

k.d  l-k,) 


If  the  retrieval  rate  decreases,  the  interval  will  become  larger  with  the  same  upper 
bound;  or,  vice  versa,  if  the  retrieval  rate  increases,  the  interval  will  become  smaller 
with  the  same  upper  bound.  However  (also  based  on  nirmerical  simulations),  there  is  a 
maximum  retrieval  rate  c* ;  if  the  retrieval  rate  exceeds  c* ,  no  control  gain  v  exists  for 
bounded  system  responses.  For  these  particular  simulation  parameters,  it  is  found 

that  c *  is  approximately  equal  to  5.l8x1d^(s ’) 


S 


’~0  OJ  02  03  0.4  05  06  0.7  OB  03  I 

/  (km) 

Fig. 3  ’’upper  a  Function  of  Tother  LengMi 

Also,  the  system  configuration  is  more  stable  with  the  long  tether  length  at  the 
initiation  phase  of  retrieval  than  witfr  the  shod  tether  length  at  the  final  phase  of 
retrieval.  Therefore,  the  nondimensional  gains,  P ,  and  p  ^ ,  are  selected  to  be  adapted 
to  the  initial  and  final  phase  of  retrieval  for  the  purpose  of  saving  control  energy,  that  is, 

g|=P,i  P2=P»i  0.1  km  -<  /  <  I  km 

(35) 

?i=g,f  km  <0.1  km 


Notice  that  in  the  case  of  retrieval  the 
nondimensional  parameters  ki .  are 
changing  with  the  instantaneous  tether 
length.  So  is  the  upper  bound  of  the  control 
gain  i>.  shown  in  Fig.  3. 


OUT-OF-PLANE  MOTION  CONTROL 

The  system  out-of-plane  motion  is  coupled  with  the  in  plane  motion  and  tether 
longitudinal  translational  motion,  especially  in  the  case  of  large  swing  angles  and 
quick  retrieval.  Therefore,  it  is  necessary  to  control  the  system  out  of  plane  motion 
during  fast  retrieval.  To  this  end,  it  is  proposed  to  control  the  shell  yaw  and  pitch  motion 
by  linear  feedback  control  law,  that  i- 

u,  =g,,vi/  +  e,,v|/•  +  g,,‘t>+g,,^• 

^  (36) 

LJ,|^  —  gjiv)/  "I"  p2»'l7  gji<t’+  Pj|*t’ 

The  selection  of  the  control  gains  p-  is  based  on  the  shell  linearized  equations  of 
motion  without  the  tethered  subsatellite  system  connected.  The  LQR  or  other  methods 
can  be  used  to  design  the  control  gains. 
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3.5  NUMERICAL  RESULTS 


First  of  all.  it  is  shown  that  with  tho  tetrieval  rate  c  =  3.45x1o'(s  ' )  (Fig. 4) 

and  without  any  controls,  the  system  responses  diverge  very  quickly  (Fig.  5). 

Fig.  6  shows  that  with  the  addition  of  the  three-axis  attitude  control  on  the  shell 
(which,  however,  does  not  involve  the  momentum  exchange  feedback  term),  the  shell 
attitude  motions  are  subjected  to  a  stable  control,  whereas  the  tether  in-plane  swing 
motion  still  diverges. 

With  the  inclusion  of  the  momentum  exchange  feedback  term  in  the  pitch  control 
law,  the  pitch  control  law  is  then  selected  as: 


Un  =  -g,e-g.0’-i)(k,a-+k,0"  ) 


.  ki-k, 

V=  — ■ — ^  -  a 
k,(l-k,) 


a  >  0 


gi=2.4  g,=  2.3  ; 

g,  =  5.4  g,=  .3.4  ; 


0. 1  km  -5  /  -5  1  km 
0.01  km^  /  <  0.1  km 


Figs.7-9  show  the  system  responses  corresponding  to  the  three  cases:  a  =  6,  a  =  4, 
and  a  =  1 .  It  can  be  seen  that  with  the  increase  of  the  control  gain  v  (a  decreases)  the 
tether  in-plane  motion  is  better  damped;  however,  the  transient  response  of  the  pitch 
motion  of  the  shell  is  subjected  to  a  larger  initial  overshoot,  but  subsequently  damps 
rapidly.  This  conclusion  agrees  with  :he  physical  explanation  of  the  momentum 
exchange  feedback  control  and  the  same  results  were  observed  in  Ref.[3).  Therefore, 
the  selection  of  the  control  gain  v  is  a  kind  of  trade-off  and  the  better  system 
responses  are  usually  achieved  with  the  control  gain  v  close  to  the  upper  bound  of  the 
stability  region. 

Fig.  1 0  shows  the  response  of  the  tether  out-of-plane  swing  motion.  The  tether 
out-of-plane  motion  will  be  disturbed  by  its  coupling  with  the  system  in-plane 
motions  and  tether  longitudinal  translation.  However,  it  can  be  seen  that  when  the 
system  in-plane  motions  are  damped,  the  tether  out-of-plane  motion  (not  directly 
controlled)  will  not  become  too  large,  with  a  maximum  swing  angle  of  20  degrees  at 
the  end  of  retrieval  (for  the  selected  control  gains  here). 

Fig.  1 1  shows  the  projection  of  the  tether  trajectory  in  the  orbit  plane. 

Fig.  12  shows  the  tether  tension  profile  during  retrieval.  It  can  be  seen  that  the 
change  of  tether  tension  is  moderate  and  no  slack  tether  occurs.  This  is  one  of  the 
advantages  of  the  current  control  method  over  tension  control. 

Fig.  13  shows  the  requirement  of  the  pitch  control  torque  during  retrieval  for  the 
selective  control  gains.  The  torque  requirement  is  reasonable. 

Although  the  control  gain  v  is  selected  in  an  adaptive  manner,  the  fluctuation  with 
the  tether  length  is  very  small  (Fig. 2).  It  will  be  very  interesting  and  of  great  practical 
value  if  a  constant  control  gain  v  can  be  applied  with  successful  results.  Fig.  1 4  shows 
some  typical  system  responses  when  a  constant  gain  v  is  used. 
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.6  CONCLUSION 


The  retrieval  dynamics  of  the  orbiting  tethered  antenna/reflector  system  without 
active  control  is  unstable.  In  the  present  paper,  the  control  concept  called  momentum 
exchange  and  compensation  has  been  applied  to  the  retrieval  of  the  tethered 
antenna/reflector  system.  The  time  rate  of  change  of  that  part  of  the  angular 
momentum  due  to  the  tether  motion  from  the  local  vertical  is  used  as  a  part  of  the 
feedback  law  in  the  antenna  attitude  control  system.  The  tether  retrieval  commanded 
length  is  taken  as  a  exponential  form.  The  selection  of  the  nondimensional  control 
gains  is  based  both  on  the  consideration  of  the  stability  of  the  system  linearized 
equations  about  a  fixed  reference  length  and  on  the  numerical  simulation  for  the  real 
time-variant  and  nonlinear  system.  The  numerical  simulation  results  have  shown  the 
effectiveness  of  the  current  control  method  for  the  suppression  of  both  the  antenna 
and  tether  in-plane  motions  at  relatively  quick  retrieval  rates. 

Another  implication  of  this  approach  to  tethered  systems  is  that  tether  swing 
librations  can  be  suppressed  through  the  control  of  the  motion  of  an  attachment  point 
or  boom  using  this  new  methodology.  This  control  concept  provides  an  alternative 
technique  for  control  of  future  possible  tethered  systems  which  may  have  some 
practical  advantages  such  as:  (1)avoiding  excitation  of  the  longitudinal  vibrations  and 
wave  propagation  of  tethers  induced  by  tension  control  and,  (2)avoiding  slack  tethers. 
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APPENDIX  A 


Q^  =  sin7  ci  -cosa  co.s7w,  +  sin7  o),.  l-sina  cosy  co. 

=  COSY  a+  cosa  sin yW:,  +  cosY  sinasin y  w, 
n,  =  Y-sina  w,-  cosa  w, 

co,,=  -  COSY  a'  -  Y'-2a(  cosasm  ycosY  (o,  +  cosYw^-sin  asm  YcosY  ) 

2  I  1  '  2  >  .-^2 

+  2Y(sinaw^  +  cosaco  J  -  ( sina  cosY  +  sinY)  co;  -cosYa)j;-(  cosacosY  +sin  y)  wj 

2 

-  2cosasin  Y  cosY  W;,  cOy-  2sinacosa  cosYw^o),  +  2sin  asin  ycosY  co^. 

2  2  ^ 

-cojj  cosa  cosY(3cos'4>  cos0-i)-6  cosa  sin  Y  cosYsin  (t>cos4)cos0 

+  6(4  sina  cosa cosYcos4)sin0cos0-  cousin Y(3sin4>cos0-l) 

+6(0^  sin  asin  YcosYsin4>  sin0cos0  -(i3^sinacosY(3  sin0-l) 

(o,,=  Y-sinao)^ -cosa ci) ,  +  s i n Y cosY a- 2 a ( cosa  cosY co ^ - s m Y cosY t^y -s m a cosY (o ^ ) 

2  2  2  .  2  .  2 
_ cosa s in  YcosYco,+sin YcosY tOy  -sinasinYcosYto,.  -  cosacos2Yto,(Oy 

'2  2  . 

+2sinacosasinYcosYa),a)j  +sinacos2Y(Oyio,  +w;  cosasinYcosY  (3cos‘(t)  cos0-l ) 

-3a)^co5acos2Ysin(t>  cos45COs"0 -(HCsina  cosa  sin  y  cosY  coset)  sin0cos0 
-coj  sin  YcosY(3sin'4)  cos0-l)+3<4Jinaco.s2Ysin<t)  sin0cos0 
+  00,.  sin  a  sm  Y  cosY  (3  sm  0  - 1 ) 

cjo  ,,  =  -cosYa  -cosasin  Ycb^-cosY(jOy  +  sinasin  Yob^+2sin  Ya  Y 

2  2  \ 

-2Y(cosa  cosYoo^-sinYojy-sinacosYwJ  +  sina  cosacosY(Wx 
-  sin  a  sin  Y  co^  cOy + cos2a  cosY  oo^  (u,  -  cosasin  Y(o,.  (o, 

_(o^  sinacosa  cosY(3cos‘(t)  cos0-l)-3odf  cos2ac(osYcos(t>sin0cos0 

_3a)^sinasinYsin<t>cos<t)cos0-3oo;cosasinYsin4)sin0cos0 

*) 

+  0)'  sina cosaco.sY(3sin0  -1) 


g^=  -sino:cosYa)^,  +  ,sin7u^+cos«  cosY ( oj‘.  +  co',  )+sin7w,io^.+  sinacosYo^^w 

1  2  2  -v  2 

4-w;  cosa  CO.S7  (3c().s4>  cos0-1) -l-3coi.. sin  7.sin4>  cos  (jocose 

-  3  w‘  s  i  11  (1  cos7  cos  ct)  s  i  II 0  coso 
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Fig.  1.  Tethered  Antenna/Reflector  System 
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Fig.  7  Tether  In-Plane  Swing  Angle  (Relative  to  Shell)  Responses 
when  Shell  Attitude  Control  is  Applied 
(Including  Momentum  Exchange  Feedback  Control  Term) 
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Fig.  8  Shell  Pitch  Angle  Responses  when  Attitude  Control  is  Applied 
(Including  Momentum  Exchange  Feedback  Control  Term) 
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Fig.  9  Tether  In-Plane  Swing  Angle  (Relative  to  Local  Vertical) 
Reponses  when  Attitude  Control  is  Applied 
(Including  Momentum  Exchange  Feedback  Control  Term) 
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Fig.  12  Tether  Tension  Profile  during  Rpirieval 


Toraue  Kml 


Fig.  13  Pitch  Torque  Requirement  during  Retrieval 
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Fig.  14  Typical  System  Responses  when  a  Cons^^nt  Gain  is  Used 


4.  THE  OPTIMAL  LQO  DIGITAL  CONTROL  OF  AM 
ORBITING  TETHERED  AMTEMMA/REFLECTOR  SYSTEM^ 


The  analysia  and  design  of  LQG  optimal  digital 
controllers  and  estimators  are  presented  for  an 
orbiting  tethered  shell  system.  Both  the  shell  and 
tether  are  considered  to  be  flexible.  In  this  paper 
the  emphasis  is  placed  on:  the  mathematical  models  of 
the  dynamics;  the  influence  of  the  number  of  actuators 
and  their  locations;  the  methods  of  measurement  of  the 
state  variables  and  design  of  LQG  optimal  digital 
controllers  and  observers. 

Attitude  and  shape  control  is  assumed  to  result  from 
point  actuators  and  tether  tension.  A  symmetric  design 
for  the  actuator  locations  is  found  by  means  of  the 
concept  of  the  degree  of  controllability  and  related 
simulations.  To  design  a  closed-loop  control  system, 
measurement  is  an  important  problem.  If  the  tether 
transverse  motion  (swing  and  transverse  vibration)  is 
not  directly  measured,  the  system  is  not  observable. 
One  of  the  main  contributions  of  this  paper  is  a  new 
practical  method  to  measure  the  tether  transverse 
motion.  With  the  measurement  of  the  tether  transverse 
motion  and  several  other  properly  placea  sensors,  the 
system  will  be  completely  observable. 

The  analysis  and  design  of  the  optimal  LQG  digital 
control  system  for  the  tethered  shell  has  been 
certified  by  simulation.  For  a  system  with  joint 
bodies  (a  flexible  shell  and  a  flexible  tether)  the 
order  of  the  system  is  very  high.  After  maximizing  the 
degrees  of  controllability  and  observability  the  best 
combination  of  the  controller  and  observer  pole 
locations  is  found  by  carefully  selecting  the 
wheighting  matrices.  Typical  figures  show  transient 
responses  under  LQG  synthesized  control  laws. 
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FIG.  1.  Tethered  Antenna/Reflector  System, 


Since  the  early  1970 's  - - -  — ell. — ==3 — ^ 

very  large  space  antennas  have 

been  proposed  for  power  ^ 

transmission,  astronomical  | 

...  rigid  boom  ; 

research  and  communications.  _[ 

The  gravity  stabilized  ,,,her  mecham.m 

configuration  is  particularly  \ 

suited  for  large  flexible  \ 

systems  to  alleviate  the  \ 

problems  associated  with  the  tether  m  > 

\  ^ 

active  control  of  large  \  \ 

structures.  A  tethered  shell  \ 

structure  (Fig.l)  to  provide  \ 

the  favorable  moment  of  \ 

inertia  distribution  for  \ 

overall  gravitational  u 

stabilization  has  been  1 

investigated  by  Liu  and  subsatTuite 

Bainum  .  In  their  paper  for  a  fig.  1.  Tethered  Antenna/Reflector  System, 

rigid  shell  and  a  flexible 

tether,  the  control  problem,  by  using  the  continuous-time 
linear  quadratic  regulator  (LQR)  technique,  has  been  solved. 
Because  of  the  huge  size  and  small  mass,  it  is  necessary  for 
the  future  analysis  to  consider  the  shell  to  be  flexible.  The 
vibrations  of  the  shell  affect  the  tether  motion  and  tension, 
and  are  coupled  to  the  rotations  of  the  shell.  To  control  the 
shape  and  attitude  of  the  shell  and  the  tether,  tether  tension 
alone  is  not  enough;  therefore,  it  is  suggested  to  add  a  few 
point  actuators.  The  number  and  location  of  the  actuators  are 
determined  by  means  of  the  concept  of  the  degree  of 
controllability.  The  state  measurement  is  another  important 
problem.  In  this  paper  it  is  shown  that  the  system  will  be 
unobservable  if  the  tether  out-of- (orbit) plane  swing  motion  or 
transverse  vibration  is  not  directly  measured.  A  method  to 
measure  the  tether  transverse  vibration  and  swing  motion  is 
developed  in  this  paper;  this  method  may  be  implemented  in 
engineering  practice,  and  can  provide  a  high  measurement 
accuracy.  It  can  be  used  for  any  orbiting  tethered  system. 
With  this  method  and  a  few  properly  placed  sensors,  it  is 


.2 


found  that  the  system  will  have  a  satisfactory  degree  of 
observability . 

In  engineering  practice  a  discrete  time  data  system  is 
more  practical  for  the  controller  based  on  the  on-board 
computer.  Thus  the  technique  of  LQG  (design  of  the  optimal 
stochastic  controller  for  a  problem  described  in  terms  of 
linear  system  models,  quadratic  cost  criteria,  and  Gaussian 
noise  models)  is  employed  to  synthesize  a  controller  with  a 
Kalman  filter  to  deal  with  the  measurement  noise  and  plant 
noise  in  the  presence  of  sampled  data  observations. 

4.2  MATHEMATICAL  MODEL 

Dynamical  Equations 

The  mathematical  model  of  a  tethered  shell  system  in  orbit 
was  developed  in  Ref.  1.  The  resulting  linearized  equations  of 
motion  were  developed  as: 

1)  shell  attitude: 

(1+QV)  <|/=  /J^c  (1) 

0''-3QV0-2Ejf"Wpn/j;-  (ii/j-;)  u'-u'po) 

+2Ei; 

JB 

d-Q’.)  ♦'-*!  1  /j; 

'  - 1  i'®,.  *  )  /  x;  “c 

2)  tether  longitude  and  swing  motions: 

<0%[rn;,(L"-u"^*3UpJ  *Ei;Ja"-3AJ  -2Ei;d„*  (2o'+20' 

-3)  i;-*m;e(3-29')  ]  {in^E,^-in,,E^)  /ntj. 

(0'43o+30)  -EH;^dJ,-3E{H;;^_^hi;j c„ 

-2E/f;^X=^ey/“c 

fll 

((|/^44)-3y)  ^ 

*3  E  (  h  i;j  b„  = 

o 

3)  tether  and  shell  vibrations: 
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+Eh;,,,(-4:-3A„)  -2ZH,j,d.*  (2a'.2e' 

-3) j;,ft  (20'-3) 

jn 

-hj;„(r+3o+30) 

E/i;^  ^  (y"-^4y-4>"-4(J))  J;„  (<|>"+4<|)-3y  )h  (gn) 

(QV3)ep_*2  Ji''”0'/M„a*  /M^a  (lon) 

0,  0:  yaw,  pitch  &  roll  angles  of  the  shell,  respectively; 
Jx,  >3y>  Jz*  inertial  moments  of  the  shell  about  the  principal 
axes  X,  y  and  z,  respectively; 

«x=  ( ny=  (•^x-‘^*)  / Jy  nz=  / ^z 

J*=Jy+m3tii2^j^  J*=j,+m3^h2+i;^ 

sip,  m^,  m^:  masses  of  the  shell,  tether  and  subsatellite, 
respectively; 

•^(  .  )  “®p-^(  .  )  /®I'  ^(  •  )  (•)  "jgZ  ^  ^  ^  ^{  .  )  (•)  “™P^(  •  )  (•)Z®Z 

=  sin(6„,/L),  ♦m=sin(nTrx/L)  the  m^^  modal  shape  functions 
of  the  tether  longitudinal  and  transverse  vibrations, 
respectively; 

ligpx,  Lgpy/  ^pz*  components  of  torque,  produced  by  Eg^  ®nd 

Ep,  which  appear  in  the  tether  force  acting  on  the  shell; 

„  f  ,  „  r  ^  [®:  tTie  external  force 

Egf.=  f  edm;  E=  edm;  E=E  ^E 

J*e  **  Jp  p  sc  acting  per  unit  mass  . 

Lgpx,  Lgpy,  Lgpxt  the  components  of  torque  contributed  by  the 

external  force  acting  on  the  shell; 
h:  the  length  of  the  boom; 
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^pn  ^n/®» 

Ap^:  the  modal  amplitude  of  the  shell; 
a:  the  base  radius  of  the  shell; 

(Oj,:  the  orbital  angular  velocity; 

the  modal  mass; 

5=r/a; 

r:  radial  distance  of  a  shell  element  from  the  symmetry  axis; 

Xpi  coordinate  of  differential  area  on  the 
surface  above  the  base  plane  of  the  shell 

*p^=ApjlaP^*CpjiP/  (RDk*pj)  ^jpikpji)  *DpjIp(Xpji)  ]cosp{&^^R^) 

the  modal  shape  function  of  shell  vibration^ 

R:  the  radius  of  the  shell  curvature; 

p,  j:  numbers  of  the  nodal  diameters  (meridians)  and  nodal 

circles,  respectively; 

Jp(.)/  Bessel  function  and  modified  Bessel  function, 

respectively,  of  the  first  kind  and  order  p; 
D=  Ep  hp3/l2(l-v2), 
bending  stiffness; 
v:  Poisson's  ratio; 

Ep;  Young's  modulus  of  the  shell; 
hpt  wall  thickness  of  the  shell; 

Wp^:  shell  frequency  of  the  n^**  mode, 

E^:  generic  force  on  the  n^**  mode; 

Fp^;  component  of  the  tether  force  acting  on  the  shell  along 
the  shell  symmetry  axis; 

♦p°:  shape  function  on  the  shell  apex; 

A„,  B„  &  C„:  the  amplitudes  of  the  longitudinal,  in-plane  and 
out-of -plane  components  of  the  m'^*'  mode  of  the  tether; 
Upg  :  the  displacement  of  shell  apex; 

<!>„;  the  n^**  mode  natural  frequency  of  tether  transverse 
vibration, 

e^^=n^n^ { l+lmg/mt) -3/4  ; 
tether  tension, 

Lgy  &  the  torques  produced  by  the  external  force; 

a,  y:  tether  in-plane  and  out-of -plane  swing  angles, 
respectively; 

Hqx»  Hgy  &  Hgx  result  from  the  external  force^; 

L  :  tether  length; 


^aui  ~ 


EAfl„(l4-m3sin2fl^/nit)  /2L 


E  :  Young's  modulus  of  the  tether; 

A  :  the  area  of  tether  cross  section; 
{  )'=d(  )/dT,  T=Ogt; 


J-2,  p«o 


P-». 


>2.  P-1 


nodal  llnaa 


It  is  noted  that  only  the 
inodes  without  nodal  diameters 
(p=0)  are  coupled  to  the 

rotation,  so  the  first  four  of  j-i,  p.o  j.2.  p.o 

them  (j“l,  2,  3  and  4)  are 

considered.  For  the  modes  with 

one  nodal  diameter  (p=l)  whose 

amplitudes  are  quite  large 

during  the  attitude  adjusting,  j-i,  p,i,  j.a 

the  first  three  of  them  (j*l,  2 
and  3)  are  included.  For  the 
modes  without  nodal  circles 

(j=0),  the  influence  of  the  .  nodai  un#* 

shell  curvature  remains  l  -  numbar  of  nodal  clrclaa 

insignificant.  Since  the  p  -  nuiUMt  ot  nodai  Mridians 

frequencies  for  these  modes  are  rig. 2  several  Modal  Shapes  of 

very  nearly  equal  to  the  a  shell 

corresponding  flat-plate 

frequencies,  the  frequencies  are  very  low.  The  mode  with  two 
nodal  diameters  is  taken  as  an  example.  For  tether  transverse 
vibration,  six  modes  are  taken,  three  for  in-plane,  three  for 
out-of-plane.  The  following  proposed  numerical  values  of  a 
tethered  shell  system  are  adopted  here: 

a  (shell  base  radius)=100  m.  lOp  (shell  mass)  =10000  kg. 

R  (shell  curvature  radius) =5  km.  H  (shell  height) =1  m. 

V  (shell  Poisson's  ratio) =1/3  hp  (shell  thickness) =1  cm. 
Ep  (shell  Young's  modulus) =8. 4047*10®  N/m^ 
m,  (subsatellite  mass) =500  kg.  m^  (tether  mass) =8. 35  kg. 
L^(commanded  tether  length) =1  km.  h  (boom  length) =80  m. 

AE  (tether  axial  stiffness) =  61645  N. 

The  deformations  of  several  modes  are  shown  in  Fig. 2,  and  the 
results  of  calculations  for  a  few  generic  modes  is  shown  in 
Table  1. 


1  -  number  of  nodal  eirelaa 

J"0,  p-2 

p  -  nuabat  of  nodal  oarldlana 
Fig. 2  Savaral  Modal  Shapes  of 
a  Shell 
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n 

(p» 

j) 

1 

“n 

^pj  ^pj 

°PJ 

T  (n) 

•^1 

1 

(0, 

1) 

3.012 

1.02778 

2.1979  0 

-.08381 

2875.55 

2 

(0, 

2) 

6.206 

1.02946 

3.1389  0 

3. 119e-3 

-322.84 

3 

(0, 

3) 

9.371 

1.03692 

3.8468  0 

-1.28e-4 

93.44 

4 

(0, 

4) 

12.53 

1.05696 

4.4425  0 

5.368e-6 

-38.73 

5 

(1, 

1) 

4.530 

1.02821 

3.8359  0 

-.01901 

0 

6 

(Ir 

2) 

7.737 

1.03198 

4.9627  0 

7.045e-4 

0 

7 

(1, 

3) 

10.91 

1.04459 

5.8752  0 

-2.845e-5 

0 

8 

(2, 

0) 

2.292 

0.00867 

3.6597  -3.47e- 

■6  0.2240 

0 

Table  1. 

The 

Result*  of  Calculation*  for  a  Few  Generic  Mode*  of 

the  Shell 

To  control  the 

Shape  and  attitude  of 

the  shell 

and  the 

tether,  tether  tension  alone  is  not  sufficient,  so  that  a  few 

point  actuators  are 

_ _ _  added.  The  concept  of 

Case  Number  Location*  of  Degree  of  degree  of  control- 

lability  has  been  used 
for  determining  the 

placement  of  the  actu¬ 

ators  on  the  shell^. 
The  degree  of  control¬ 
lability  is  the  scalar 
measure  of  system  con¬ 
trollability  and  its 

reciprocal  indicates 
the  effort  to  control 
the  system.  The  values 
of  the  degree  of  con¬ 
trollability  for  some 

actuator  placement 
designs  are  listed  in 
Table  2.  It  is  assumed 
that  the  thrusters  at 
the  shell's  edge  have 

two  jet  directions: 

Table  2.  The  Degree  of  Controllability  for  a  i)  tangent  to  the  edge; 

Few  Actuator  Deaign*  ii)  normal  to  the  shell 

surface.  Each  of  the 
other  thrusters  has 


2 

3 


6 


12 

12 

12 

12 

8 


0.0 

1.457*10‘''’ 

0.0 

0.0 

5.893*10'” 

5.443*10*’ 

5.443*10'’ 

4.987*10’’ 
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only  one  jet  direction,  i.e.  normal  to  the  shell  surface.  In 
case  8  only  eight  actuators  are  used,  and  the  degree  of 
controllability  is  satisfactory,  so  this  actuator  location 
design  is  recommended. 


For  an  actuator  which  can  generate  a  force  fy  f^)'^ 
and  is  placed  at  a  location  y^  control  torque 
is  given  by  x  so  in  eqs.  (l)-(3} 


Lepy  X 


^py  ^pxJ  ^d®pz  ^d®*py)  (12) 

where  is  the  distance  between  the  mass 
centers  of  the  shell  and  the  tethered 
shell  system.  In  the  remaining  equations 

L^y  “  -Jjf  Epg  /  /Hj.  (13) 

<151 

(16) 

(17) 

E„-t  ♦,.(y,.z,)f,  (18) 

(each  jet  is  along  Che  shell's  x  axis) 


1/  double 
^  inner 
frame 

A, -A. 

bearings  & 
angular  tether 
transducers 


The  system  will  be  unobservable  if 
tether  transverse  motion  (swing  and 

vibration)  is  not  measured.  A  new  method  subsatellite 

to  measure  the  motion  is  shown  in  Fig.  3a.  , 

Fig. 3a.  A  Design  for 

The  accuracy  of  angular  transducer  can  be  Measurement  of 

very  high  (e.g.  the  ones  used  on  inertial  of*"the'^*Tether'’" 

platforms)  during  station  keeping  (i.e. 
the  error  angle  is  very  small).  The 

subsatellite  is  quite  heavy,  so  that  an  inner  framed  structure 
is  adopted.  Since  the  tether  can  bear  only  the  longitudinal 
tension,  the  position  of  the  frame  should  coincide  with  the 
tangent  direction  of  the  tether  at  the  attachment  point 


(Fig. 3b).  The  displacement  of  the  tether  in-plane  vibration  is 
approximated  by 

m 

W=T.C^{t)  sinCnnx/L)  then 

n*X 

o.-tana.=df//dxL.o*EC-(c)  (nn/Dcosinnx/ L)  |^o=E  C„{t)  inn/L) 

n«l  /i«l 

The  angular  transducer  A2  is  assumed  to  measure  the  in-plane 
displacement  between  the  shell  reference  axis  and  the  tether 
line 

yi=Ar,  (oe+ffi)  [«+  E  C„(  t)  (nn/L)  ]  (19) 

n«l 

Similarly  for  out-of-plane  case,  the 
value  obtained  from  transducer  is 

(nit/L)]  (20) 

A  dial  with  magnetic  encoder  is 
assumed  to  be  used  to  measure  tether 
length  and  its  changing  rate.  The  dial 
is  fixed  to  the  reel  from  which  the 
tether  is  retrieved  or  deployed.  When 
the  dial  turns  it  produce  a  series  of 
impulses.  The  sum  of  the  impulse  number 
represents  tether  length,  and  the 
reciprocal  of  impulse  interval 
represents  the  rate  of  tether  change. 

y3=/C,€  (21)  74=^  (22) 

where  e=(L-Lg)/Lg 
L,;  is  the  commanded  length  of  the  tether. 

The  shell  attitude  is  assumed  to  be  measured  by  an 
inertial  platform  based  on  gyroscopes.  The  advantage  of  the 
inertial  platform  is  that  it  produces  signals,  for  both  the 
attitude  and  angular  velocities,  with  a  high  signal: noise 
ratio.  The  drift  of  the  inertial  devices  is  then  modified  by 
two  infra-red  horizon  sensors  of  the  Earth  and  two  sun 
sensors^.  The  reason  why  an  inertial  platform  is  used  is  that 


Fi9.3b.  The  In-plan« 
Angl*  Between  the 
Deformed  Tether  and 
the  Shell  Reference 
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the  infra-red  and  sun  sensors  give  signals  with  low  ratios  of 
signal  to  noise,  and  in  the  dark  part  of  the  orbit  the  sun 
sensor  can  not  give  any  signal  at  all.  So  the  measurement 
equations  for  the  shell  attitude  are: 

(23)  ♦  (24)  y-,=K^^  (25) 

(26)  (27)  (28) 

It  is  assumed  that  four  displacement  sensors,  which  are 
collocated  with  the  actuators  in  the  first  and  second 
quadrants,  are  used  to  measure  the  shell  vibration: 

(5i  Ci)ep^(t)  (i=ll,  12, 13, 14)  (29) 

Eqs. (19) -(29)  consist  of  the  measurement  equations,  for 
which  the  system  is  completely  observable. 

The  State  Bquatloas  and  LQd  Regulators  and  Observers 

After  substituting  Eqs.  (ll)-(15n)  into  Eqs.  (l)-(lOn), 
and  after  some  algebraic  manipulations,  the  system  equations 
can  be  written  in  the  state  vector  form: 

x'  »  A  X  +  B  u  (30) 

where  u-  (f^,  fj, . fg,  Af)  ^ 

(®'  . <73,  Y,i3i,b2,i?3, 

^ ,  a',e',ej^,ej^, . ,€p^,  c[,  Cj,  c,',  ,y' ,b[,b2.b!i)  ^ 

Ci^CjLc  b^^BjL^  Af»F„/  ((i)\fl7;cLp)  +3  [ 

In  fact,  dynamic  systems  are  driven  not  only  by  their  own 
control  input,  but  also  by  disturbances  which  we  can  neither 
control  nor  model  deterministically.  Sensors  generally  do  not 
provide  exact  readings  of  desired  quantities,  but  introduce 
their  own  system  dynamics  and  distortions  as  well.  Furthermore 
these  devices  are  also  noise  corrupted.  In  order  to  solve 
these  difficulties,  the  LQG  techniques  (design  of  the  optimal 
stochastic  controller  for  a  problem  described  in  terms  of 
linear  system  models,  quadratic  cost  criteria,  and  Gaussian 
noise  models)  are  effective.  In  this  section  the  LQG  theory  is 
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applied  to  the  design  of  suboptimal  regulators  and  estimators 
for  the  tethered  shell  system. 

For  the  infinite-time  optimal  regulator,  given  the  system 
represented  by  Eqs.(30),  the  task  is  to  find  the  optimal 
control  such  that  the  performance  index 

J^\r{x^O^*u'^R^u)dt  (31) 

A  J  O 

is  minimized.  In  Eq.(31)  is  a  symmetric  positive 

semidefinite  matrix;  is  a  syTnmetric  positive  definite 

matrix.  After  discretization  and  consideration  of  noise. 
Eq.(30)  becomes^ 

x(k+l)  =A(T)x{k)  *B(T)  u{k)  *w(k)  (32) 

y  (k)  =  C  X  (k)  +  V  (k)  (33) 

where  ^CD-e*^  B(T) -f  ^e*'^  B  dt 

J  O 

and  Eq.(33)  is  the  state  vector  form  of  Eqs. (19)-(28) . 
w(k)  and  v(k)  are  the  n-dimension  and  m-dimensional  white 
Gaussian  discrete  time  noise  with 

E[w{i)  w*’(j)  1  =06^^  Bivii)  v^{j)]-Rh 

respectively,  and  assumed  to  be  independent  of  each  other  and 
the  initial  conditions.  The  initial  condition  x(0)  is  modelled 
as  a  Gaussian  random  vector  with  mean,  x^  and  covariance,  P^. 

From  Eq.(31),  the  discretized  performance  index  is 

j-.Af  [x^{k)Cx(k)^2x'^(k)ylu(k)^u^(k)^u(k)]  (34) 

2  jt.o 

where  (?  *  J  0^  e**^  dt  -  J  ^  e*^^  gi  t,  0)  dc 

R^jJ[Rc*9^<t:>0)Ogg{t,0)dC  g(t,0)  B  ds 


The  problem  now  is  to  find  the  control  u(k)  ,  so  that  the 
performance  index  in  Eq. (34)  is  minimized.  If  the  system 
represented  by  Eq. (32)  is  controllable  (or  stabilizable)  and 


observable  (detectable),  the  LQG  solution  is^'® 

u(k)  *  -G  Ji(k)  (33) 

G=  (34) 


where  P  satisfies  the  algebraic  Ricatti  equation 

P=A'^PM0-  (A'^Pl^+f^  (if+^^PS) -MJ’^PA+W'^)  (35) 

The  state  estimate  at  the  interval  can  be  related  to  its 
predicted  value  based  on  information  from  the  (k-1)*^  interval; 
together  with  the  measurement  vector  at  the  k^^  interval,  as 

JUk)  »J?(Jc/ic-l)  *K[y{k)  -C^ik/k-1)  ]  (36) 

where  ^(k/k~l)  •A  ^(k-1)  *  B  u(k-l)  (37) 

k-  *  P.  C’’  (  C  P.  (:*■  +  P  )-i  (38) 


The  covariance  of  the  state  estimate,  P^,  satisfies  the 
algebraic  Ricatti  equation 

P^^{A-K’C)  P^(A~K’a  ^*K*RK''^-^0  {  where  K*^AK  ) 


Due  to  computational  delay,  what  is  often  considered  for 
implementation  is  not  the  control  law  given  by  Eq.(33),  but  a 
subopt imal  control 

u(k)  =  -GP(k/k-l)  where 

(40) 

k(k*l/k)  =JUt(k/k-l)  *Su{k)  +k’[y(k)  -CP(k/k-l)  ] 


It  is  wise  to  seek  the  possibility  of  a  separate  determination 
of  the  state  estimate  and  the  controller.  The  error  of  the 
predicted  estimate  of  the  state  is  defined  as 

Ak(k/k-l)=r(k)-P(k/k-l)  (41) 

the  x(k*»>l)  and  Ax(k+l/k)  may  be  combined  into  a  simple  system 
described  by  the  state  equations^ 


x(k*l) 

A~Sg  BG 

x(k) 

J  0 

^(k) 

d3if(k+l)  /k) 

0  A-K*C_ 

Ax(k/k-l). 

▼ 

-T 

y(k) 

From  Eq.(42),  the  LQG  suboptimal  control  dynamics  consist  of 


two  parts;  the  dynamics  of  the  plant  with  a  feedback 
controller  and  the  dynamics  of  the  estimator  feedback  loop. 
The  matrices  describing  the  dynamics  of  the  closed  loop 
controller  and  the  closed  loop  estimator  are 

A*=A-BG  (43)  and  A*=A-X*C  (44) 


respectively.  The  eigenvalues  of  A^*  and  A,*  are  independent  of 
each  other.  In  the  synthesis  it  is,  therefore,  possible  to  ar¬ 
range  the  poles  of  the  estimators  and  controllers  separately. 

It  is  well  known  that  if  and  only  if  the  open  loop  system 
model  is  both  stabilizable  with  respect  to  u(k)  and  detectable 
with  respect  to  y(k),  there  do  exist  gains  G  and  K*  that  can 
provide  asymptotic  closed  loop  stability.  Furthermore,  under 
the  stronger  assumptions  of  complete  controllability  and 
complete  observability  (as  in  the  designs  of  the  actuators  and 
measurements  for  the  tethered  shell  system,  here) ,  we  can 
place  both  the  regulator  and  observer  poles  arbitrarily 
(within  the  restriction  of  the  complex  conjugate  pairs) . 

It  should  be  noted  that  in  Eq.(30)  the  differentials  are 
respect  to  the  orbit  time,  i.e. 


'»dXi/d«gt=  ( l/«g)  dx^/dt 


If  the  noise  added  to  dx^/dt  is  w^,  the  noise  added  to  x^^' 
should  be  (1/Uq)Wj^.  So  the  covariance  matrixes  should  be^°: 


0 

(l/«’c)  J20 


0 


(45) 


For  the  control  loop  let  the  weighting  matrices 
Qc-C^  C 

In  order  to  find  the  appropriate  arrangement  of  the  observer 
and  controller  poles  it  is  necessary  to  study  the  loci  of  the 
eigenvalues  of  A,*  and  A^.*  with  Mq,  and  respectively.  The 
maximum  and  minimum  moduli  of  the  eigenvalues  of  A,*  and  A^* 
vs.  the  different  parameters  ^g,  and  are  listed  in  Table 


3  and  4,  respectively. 

Based  on  these  data  in  Tables  3  and  4  we  can  arrange  the 
position  of  the  controller  and 

observer  poles.  Table  3.  Max.  &  Min.  Moduli 

of  Eigenvalues  of  A  *  vs. 

Mg  (Aig»10'®) 

4.3  SINULATIOMS 


It  is  assumed  that  the 
accuracy  of  the  displacement 
sensors  is  about  1  cm.  for  the 
shell  deflection,  the  angular 
sensor  accuracy  is  about  10~^ 
rad. ,  and  the  modelling  error  for 
the  dynamical  system  is  less  than 
the  error  of  the  measurement 
sensors.  In  Eq. (45)  the  parameters 
of  the  covariance  matrices  used  for 
simulation  are  =10"®  and 
It  is  also  assumed  that  both  the 
measurement  and  the  plant  noises  are 
Gaussian  white  noises  with  zero 
mean. 


Pg  min.  max. 

10'”  6.35140-4  0.96714 

10‘’2  1.68170-2  0.97003 

10"’^  0.11649  0.97714 

10"’®  0.55691  0.98687 


Tablo  4.  Max.  &  Min Moduli 
of  Eigonvaluos  of  A^*  vs. 


min.  max. 

10*®  0.88481  0.99948 

10*‘  0.83930  0.99943 

10*®  0.73230  0.99850 

10*'  0.57579  0.99701 


In  general,  as  for  the  design  of 
the  Kalman  filter,  the  considered  covariances  of  the  measure¬ 
ment  and  plant  noises  should  be  greater  than  or  equal  to  the 
actual  ones.  After  these  conditions  have  been  satisfied  the 
location  of  the  observer  poles  may  be  changed  by  means  of  the 
variation  of  the  parameters  and  Hq,  When  4j^=*l0"®,  4^=10"^® 
the  minimum  and  maximum  moduli  of  the  eigenvalues  of  A^* 
(observer)  are  0.55691  and  0.98687. 

As  we  know  the  minimum  and  maximum  moduli  of  the  observer 
closed-loop  eigenvalues  should  be  less  than  those  of  the 
controller  closed-loop  eigenvalues,  so  that  the  estimator  can 
provide  accurate  timely  state  information  for  the  controller. 
Based  on  this  principle,  the  possible  appropriate  minimum  and 
maximum  controller  eigenvalues  are  0.57579  and  0.99701. 

The  initial  state  is  also  assumed  to  be  noise-corrupted 

x(0)=*o+Wo  Xo«E{*(0)}  Po=E{Wo  w/} 
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The  initial  conditions  are  assumed  as: 

(0) =0 (0) =0 . 1  rad.  Xj  (0)  =a (0) =0 . 1  rad  Xj  (0) =e (0) =0 . OOi 
x<  (0)  =e  (0)  =0 . 001  Xij{0)  =0^(0)  =0.001  Xi5  (0)  =i|;  (0)  =o.  1  rad. 
x^g  (0)  =<^(0)  =0 . 1  rad.  Xi7  (0)  =Y  (0)  =0. 1  rad.  x^g  (0)  =hi (0)  =0  .  OOi 

and  the  rest  of  the  initial  conditions  are  taken  to  be  zero. 
The  initial  values  of  the  estimated  state  are  assumed  to  be 
zero,  i.e., 

Xi^(0/-1)=0.0  (i=l,  2,  . .  40) 

The  strategy  of  choosing  the  sampling  time,  T,  should  be 
to  select  the  sampling  time  as  long  as  possible  after  the 
performance  of  the  sampled-data  system  meets  the  requirements 
of  the  design.  Generally  speaking,  too  long  a  sampling 
interval  tends  to  deteriorate  the  performance  of  a  sampled 
data  system  (increase  the  sensitivity,  decrease  stability, 
loss  of  controllability  etc.)^°  On  the  other  hand,  the 
implementation  of  a  very  short  sampling  interval  may  be 
limited  by  computer  operation  times  and  the  expense  of  fast 
A/0  and  D/A  devices.  The  choice  of  the  sampling  time,  T,  will 
also  be  constrained  by  the  capability  to  restore  the  signal 
which  comes  from  observational  data,  i.e.  T  meets  Shannon's 
theorem,  t<Tr/(^g,  where  is  the  maximum  frequency  of  the  input 
signal.  In  the  simulation  here  T-5  sec. 

The  result  of  the  simulations  are  shown  in  Figs. 4-6.  Fig. 4 
is  for  the  in-plane  motion.  Fig. 5  for  the  out-of -plane  motion, 
and  Fig. 6  for  each  actuator  force  needed  during  the  regulating 
process.  Since  the  tether  can  not  support  compressive  forces 
and  a  large  tension  force  should  be  avoided,  the  tether 
tension  boundary  is  set;  0<Ftx<  5  newton  (Fig. 4  c)  .  The  result 
of  the  simulation  is  satisfactory.  The  initial  errors  are 
converging  smoothly,  the  required  forces  for  each  actuator  are 
reasonable.  After  5  minutes  we  can  hardly  see  any  difference 
between  the  actual  value  of  each  state  variable  and  its 
estimate. 
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CONCLUSIONS 


1)  The  actuator  placement  design  used  here  is  satisfactory. 
The  degree  of  controllability  is  quite  high  considering 
only  eight  point  actuators  are  applied.  This  design  can 
control  elastic  vibrations  and  orientations  of  a  flexible 
shell  and  a  flexible  tether,  and  it  has  been  verified  by 
simulations. 

2)  The  system  is  not  observable  if  the  tether  transverse 
motions  are  not  directly  measured.  A  new  method  of 
measuring  tether  transverse  motion  presented  here  can  be 
applied  not  only  to  tethered  shell  systems,  but  also  to  any 
system  with  a  tether.  This  design  is  easy  to  implement  in 
engineering  practice  and  of  high  accuracy  of  measurement. 

3)  The  discrete  time  data  system  is  more  practical  for  a 
controller  based  on  an  on-board  computer  system.  The 
sampling  time,  T»5  sec.  is  a  good  compromise  between  the 
performance  of  the  system  and  the  capacity  of  the  on-board 
computer . 

4)  For  a  dynamic  system  with  plant  and  measurement  noises,  the 

LQG  technique  is  effective.  As  for  the  placement  of  the 
controller  and  observer  poles,  the  minimum  and  the  maximum 
moduli  of  the  eigenvalues  of  the  closed-loop  observer  must 
be  less  than  the  minimum  and  maximum  moduli  of  the 
eigenvalues  of  the  closed-loop  controller,  respectively,  so 
that  the  observer  can  provide  the  timely  and  accurate 
estimate  of  the  state  variables  for  the  controller.  We 
should  also  ensure  that  the  ratio  can  not  be  too 

small,  otherwise  the  Kalman  filter  may  become  too  sensitive 
to  the  observational  noise. 

5)  Since  the  tethered  shell  system  consists  of  a  flexible 
shell  and  flexible  tether,  the  dimensionality  of  the  state 
vector  is  as  high  as  40.  Because  of  the  practical  pos¬ 
sibility  of  on-board  computational  implementation,  it  is 
suggested  for  further  work  to  design  low-order  controllers 
based  on  robustness  theory. 
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Fig. 4. The  transient  response  of  the  in-plane 
motion  of  a  tethered  shell  system 
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5.  EFFECTS  OF  SOLAR  RADIATION  PRESSURE  ON  THE  TETHERED 
ANTENNA/ REFLECTOR  SUBSATELLITE  SYSTEM 


5.1  Introduction 

since  the  initial  proposal  of  the  tethered  subsatellite  system, 
the  dynamics  and  control  of  the  tethered  subsatellite  systems  have 
been  investigated  by  a  host  of  investigators.  A  comprehensive 
survey  article  was  given  by  Misra  and  Modi  recently.  Sources  of 
applications  of  tethers  in  space  have  been  proposed  for  different 
purposes.  For  example,  a  number  of  very  large  space  antenna  and 
reflector  types  of  orbiting  structures  have  been  proposed  for  power 
transmission,  astronautical  research  and  communications. 

A  tether  tension  control  law  based  on  the  tether  length  and 
length  rate  for  in-plane  control  was  formulated  by  Rupp  f  ^ 1  and 
represented  a  pioneering  effort  in  establishing  the  feasibility  of 
control  by  using  tether  tension  control  modulation.  Recently, 
Bainua  and  Liu  developed  a  system  mathematical  model  of  a  class 
of  large  antenna  and  reflector  orbiting  structures  which  has  an 
articulated  tether  connected  supporting  mechanism  to  provide  the 
favorable  moment  of  inertia  distribution  for  over-all  gravitational 
stabilization  and  controllable  tether  tension  used  for  producing 
restoring  torques  on  the  shallow  shell  reflector.  The  optimal 
control  law  for  tether  tension  control  has  been  suggested  and 
evaluated.  The  numerical  results  show  that  the  in-plane  motion  of 
the  tether  could  be  asymptotically  stable  with  Rupp's  tension  law; 
and  the  transient  response  can  be  improved  significantly  by 
carefully  selecting  optimal  control  gains  with  suitable  state 
penalty  and  control  penalty  matrices. 

As  the  different  control  laws  are  being  developed,  a  question 
arises  that  if  some  environmental  disturbances  are  considered,  will 
the  tension  control  laws  still  be  suitable  for  maintaining  the 
required  pointing  accuracy  of  the  antenna /ref lector  subsatellite 
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system  ? 

•  It  is  known  that  the  solar  radiation  pressure  is  one  of  the 
dominant  environmental  disturbances,  especially  in  higher  orbits. 
An  overview  about  the  environmental  torques  was  given  by 
Shrivastava  and  Modl.^^^  The  tension  control  is  based  on  utilizing 

•  the  gravity-gradient  and  centrifugal  forces;  therefore,  the 
capability  of  the  tension  control  will  be  decreased  as  the  altitude 
level  is  increased.  On  the  other  hand,  the  solar  radiation 
disturbance  will  be  augmented  because  of  less  oxygen,  nitrogen  and 

•  ozone  at  the  higher  altitude. 

In  general,  the  disturbance  forces  and  torques  induced  by  solar 
radiation  pressure  vary  with  the  illuminated  area,  the 
characteristics  of  the  illximinated  surface  of  the  orbiting 

•  structure,  and  the  offset  of  the  center  of  solar  radiation  pressure 
from  the  center  of  mass  of  the  system.  The  general  formulations  for 
evaluating  the  solar  radiation  pressure  forces  and  torques  are 

as  follows: 


■ 

(5.1) 

• 

Fj^=-2wf  {]fi-d)^]^d3 

(5.2) 

£=  (1-e) 

(5.3) 

• 

(/i-d)zxdds 

(5.4) 

• 

N^=-2  wj  ( i?-d )  ^zx/ids 

(5.5) 

(l-e)^+€^ 

(5.6) 
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where 


resultant  solar  radiation  pressure  force  on  a 
completely  absorbing  surface 
ZfS  resultant  solar  radiation  pressure  force  on  a 

completely  reflecting  surface 
F:  resultant  solar  radiation  pressure  force  on  a 

general  surface  with  arbitrary  reflectivity 
Ha 5  resultant  solar  radiation  pressure  torque  on  a 

completely  absorbing  surface 
Hf:  resultant  solar  radiation  pressure  torque  on  a 

completely  reflecting  surface 
H:  resultant  solar  radiation  pressure  torque  on  a 

general  surface  with  arbitrary  reflectivity 
e;  reflection  coefficient  i.e.  the  ratio  of  the  energy 
density  of  the  reflected  light  flux  to  the  energy 
density  of  the  incident  light  flux 
s:  illuminated  part  of  the  surface 

But  there  is  a  major  difficulty  in  determining  the  boundary  of 
the  illuminated  area  on  the  surface  of  the  orbiting  structure.  The 
illuminated  area  is  affected  by  the  position  of  the  orbiting 
structure  in  the  orbit,  the  attitude  and  the  shape  of  the 
structure.  Sometimes  it  is  extremely  difficult  to  obtain  an  exact 
analytical  formulation  for  the  solar  radiation  disturbance.  Bainum 
and  Krishna  used  a  numerical  method  to  approximate  the  force  and 
moment  induced  by  the  solar  radiation  pressure  on  an  orbiting 
flexible  beam  and  plate,  because  of  the  complicated  vibrational 
modes  and  the  time  varying  shadow  attributed  to  each  of  those 
modes . 

One  of  the  objectives  of  the  study  is  to  develop  the  analytical 
formulations  for  the  modelling  of  the  solar  radiation  disturbance 
on  a  class  of  large  tethered  antenna /ref lector  types  of  orbiting 
structures.  The  other  is  to  evaluate  the  effects  of  these 
disturbances  in  the  case  of  station-keeping  and  to  determine  at 
which  altitude  level  the  tension  control  law  is  still  able  to 
maintain  tb  >  satisfactory  pointing  accuracy  of  the  tethered 
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antenna/reflector  subsatellite  system  by  numerical  simulation. 

5.2  Modelling  of  the  Solar  Radiation  Disturbance  Torques 

The  system  discussed  here  is  composed  of  a  spherical  shallow 
shell  reflector,  a  boom  attached  at  the  apex  of  the  shell 
reflector,  an  articulated  tether  connected  at  the  end  of  the  boom 
and  a  spherical  subsatellite. (Fig.  5.1)  For  mathematical  modelling 
of  the  solar  radiation  disturbance,  the  following  assumptions  are 
made: 

1)  The  shell  reflector,  boom  and  the  subsatellite  are  assumed  to 
be  rigid  bodies,  and  the  tether  is  compressible  but  considered 
as  a  straight  line  (neglecting  the  transverse  deformation) 
during  the  swing  motion.  The  surfaces  of  these  bodies  are 
characterized  by  isotropical  reflectivity  properties. 

2)  The  spherical  shallow  shell  is  totally  illuminated  either  on 
the  inner  shell  surface  or  outer  shell  surface,  since  the 
ratio  of  the  height  of  the  shell  to  the  radius  of  the  base  of 
the  shell  is  far  less  than  unity. 

3)  Local  shadowing  on  the  tether  and  the  boom  is  considered. 
There  are  two  phenomena;  one  is  that  some  sections  of  the 
tether  or  the  boom  will  not  be  illuminated  in  certain 
positions  of  the  orbit  due  to  the  shadowing  of  the  reflector; 
the  other  is  that  the  boom  or  the  tether  is  considered  as  a 
cylindrical  column;  therefore,  at  any  instant  of  time  the 
tether  or  the  boom  can  only  be  illuminated  on  one  side  when 
the  emission  of  the  sun  light  and  the  reflected  sun  light  from 
the  Earth  or  other  space  objects  are  neglected.  In  addition, 
the  shadow  of  the  Earth  on  the  orbit  has  been  taken  into 
account  during  the  numerical  simulation. 

4)  The  solar  incidence  vector  in  the  orbital  coordinate  reference 
system  is  treated  as  a  time  varying  parameter.  It  varies  with 
the  inclination  of  the  orbit  and  the  position  of  the  space 
structure  in  the  orbit.  A  circular  and  equatorial  orbit  is 
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assumed  in  the  numerical  simulation  without  losing 
generality. 

The  coordinate  systems  used  in  the  development  of  the  system 
equations  of  motion  are  shown  in  Fig.  5.2;  OXqYqZq  is  an  orbital 
reference  system  centered  at  the  center  of  mass  of  the  shell 
reflector  with  OX^  along  the  local  vertical,  OY^  along  the  normal 
to  the  orbit  plane  and  OZ^  along  the  orbital  tangent  velocity 
direction.  OXj.Yj.Zp  is  a  shell  body  reference  frame,  the  axes  OXp, 
OYp  and  OZp  are  principal  axes  of  the  shell  reflector.  O^X^Y^Z^  is 
the  tether  reference  frame  with  O^X^  along  the  undeformed  tether, 
where  0,.  is  the  point  from  which  the  tether  is  deployed  or 
retrieved.  The  coordinates  of  0^  in  the  shell  reference  are 
(h^,0,0)  . 

The  yaw,  pitch  and  roll  angles  of  the  shell  are  denoted  by  i|t,  6 
and  <p,  respectively.  An  Euler  angle  rotation  sequence  is  assumed 
as;  (l)t,  (2)0  and  (3)0.  Therefore  the  transformation  from  OX^Y^Zq 
to  OXpYpZp  is  assumed  to  be  given  by; 


where 

O0C0  54>CTjr+o0s0si)r  s0s0-oi|>s0CT|r 
=  -s^c6  ci(>a|r-s4>s0si|r  oi^siir+5<t>s6ci|r  (5.8) 

s0  -c6si(r  c6c!|r 

The  in-plane  swing  angle  of  the  tether,  a,  depicts  the  angle 
between  the  projection  of  the  non-deformed  tether  in  the  orbit 
plane  and  the  OXp  axis;  and  the  out-of-plane  swing  angle  of  the 
tether,  y,  indicates  the  angle  between  the  projection  of  the 
undeformed  tether  in  the  OX^Yp  plane  and  the  OXp  axis.  The 
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transformation  from  O^X^Y^Zt  to  OX^Y^Zj.  can  be  derived  by; 


rt 

= 

0 

Zr. 

0  . 

(5.9) 


where 


Ot  = 


-cyca  sy  sacy 
syca  cy  -sasy 
-sa  0  -ca. 


(5.10) 


The  unit  solar  incidence  vector  can  be  formed  in  the  orbital 
reference  system: 


d,  =  [a,,  b,.  c^y 


(5.11) 


where  a^,  b^  and  c^  are  components  along  the  OXq,  OYq  and  OZq  axes, 
respectively.  Furthermore,  the  unit  solar  incidence  vector  can  be 
#  transformed  into  the  shell  reference  system: 


d 


(5.12) 


here  a^,  b^  and  c^  represent  the  components  along  the  OX^,  OY^  and 
OZ^  axes,  respectively.  Similarly,  the  unit  solar  incidence  vector 
in  the  tether  reference  system  can  be  written  as: 


d 


(5.13) 
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a^,  and  c,.  are  the  components  of  the  unit  solar  incidence 
vector  along  O^X^,  and  O^Z^.  axes,  respectively. 

The  disturbance  torque  induced  bv  the  solar  radiation  pressure  on 
the  antenna /reflector  can  be  expressed; 

1)  for  completely  absorbing  shell  surface 


0  ■ 

= 

= 

2) 


for  completely  reflecting  shell  surface 


'  0 

^rfy 

^rtz. 

»  R, 

<^r 

(5.14) 


(5.15) 


3) 


for  general  shell  surface  with  arbitrary  reflectivity  properties 


where 


= 


N, 


ry 


=  {l-e)N^*eN^  , 


J?,=  -nwR^c,sin2Po  (I-cosPq--^) 

I\ 


R,=  -xwi?^c,sin*Po  (1--^) 

A 


(5.16) 
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w;  solar  constant 

Sq:  cone  angle  of  the  spherical  shell  reflector 

R:  radius  of  the  spherical  shell  reflector 

rQ:  distance  from  the  apex  of  the  shell  to  the  center  of 

mass  of  the  shell 

The  disturbance  toroue  induced  bv  the  solar  radiation  pressure  on 
the  boom  has  the  form: 

1)  for  completely  absorbing  boom  surface 


0 

^  = 

II 

2)  for  completely  reflecting  boom  surface 


0 

^  = 

Ntfy 

=  Be 

--^sinUi,--|-b,c,cos^l^-c|  (sinA.^-  ) 

_^bfz 

( cos  A,  + -|brC,sinU  cosU  ^ 

(5.18) 


3)  for  general  boom  surface  with  arbitrary  reflectivity  properties 


where 


B,=-wr*(h,cosX^+c^sinX^)  dl-llo) 


Bf=-2wr^(ll-llo) 


r^:  radius  of  the  cylindrical  column  boom 

1|,:  length  of  the  boom 

I^q:  length  of  the  shadowed  boom 


^  b 


tan'^-^ 

(VO) 

tan'^-r^  +  11 

(VO) 

The  disturbance  force  induced  bv  the  solar  radiation  pressure  on 
the  subsatellite  is  written; 

1)  for  both  absorbing  and  reflecting  subsatellite  surface 


F 

at 

^afx 

F 

*  smy 

*  nwrj 

i>t 

=£^  = 

^aty 

F 

aazj 

,^sfz. 

2)  for  general  subsatellite  surface  with  arbitrary  reflectivity 
properties 

=  (1-e)  F^.*eF,r  ,  (5-21) 


here  rg  is  the  radius  of  the  spherical  subsatellite. 


the  subsatellite  about  the  point  is  written: 


1)  for  both  completely  absorbing  and  reflecting  subsatellite 
surface 


ij-wnrl  ~^t 


sfyi  / 


(5.22) 


2)  for  general  subsatellite  surface  with  arbitrary  reflectivity 
properties 


where  1^:  the  instantaneous  length  of  the  tether 


(5.23) 


ure  on 


1)  for  completely  absorbing  tether  surface 


=2wr^(b^C03X^*C^3ink^)  .  (5-24) 


2)  for  completely  reflecting  tether  surface 


^ tf~  ^cfy  ~  4  Wit  (  i  g-1  j(j) 


bl  (cosX  c-  ^ )  +-^cos^X  j+ybjCcSin^A.  ^ 
-^sin^X  (sinX  )  +-|-jbj.Cj.cos^X  ^ 


(5.25) 
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3)  for  general  tether 
properties 


surface  with  arbitrary  reflectivity 


tx 


ty 


tt 


(5.26) 


where  r^:  the  radius  of  the  tether 

l^g:  the  length  of  the  shadowed  tether 


tan*^-^ 

(b,>0) 

tan'^~  +  n 

(b,<0) 

The  disturbance  torque  induced  bv  solar  radiation  pressure  on  the 
tether  about  the  point  0^  is  given  by; 

1)  for  completely  absorbing  tether  surface 


■  0 

^  = 

=2wrc(27j.cosX  c+CpSinX  j)  ill-llo) 

2)  for  completely  reflecting  tether  surface 


tfx 


tty 


tfz 


=  2wz 


sinU, 


■  sin^A.  J.-CI  (sinX 


3  )  --h,c,cos^X  , 


bl  (cosX  c-  ^  )  +-^cosU  £+-^bj,C(.sin^X 


3)  for  general  tether  surface  with  arbitrary  reflectivity 
properties 


N.  = 


N, 


tx 


ty 


tz 


(5.29) 


5.3  Modelling  of  the  System  Equations  of  Motion 

since  the  above  analytical  models  of  the  solar  radiation 
disturbance  are  derived  without  taking  the  vibrations  of  the  tether 
and  the  shell  into  consideration,  the  system  equations  of  motion 
must  also  be  developed  based  on  this  assumption. 

Besides,  it  is  known  that  the  roll  and  yaw  motion  of  the  shell  as 
well  as  the  out-of-plane  swing  motion  of  the  tether  are  decoupled 
from  the  pitch  motion  of  the  shell  and  the  in-plane  swing  motion  of 
the  tether  in  the  linear  range.  (^1  Thus,  after  placing  the  solar 
radiation  disturbances  derived  previously  into  the  right  hand  side 
of  the  equation  of  motion  and  proceeding  with  a  series  of 
complicated  algebraic  manipulations,  the  linear  non-dimensional 
equations  for  in-plane  motion  are  obtained: 


U>Uy 


(5.30) 


+8'' + 3  0 + 3  a -ic4e^c= 


(5.31) 


ec+2  (k+P  je'-3ke,>2ka'=3  (k+p  J  + 


F  F 

^  sx  ^  cx 


SC 


(5.32) 
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where 


m. 


m^f^—iTig+nig.  i  ic=- 


m 


St 


kj_=- 


2IJtc 


^2=-Jx-7  '  ^4=  rj  >  Jc5=-2/n,^^-^ 


j: 


Jy  =  Jy*hxnigt;~^X^X  •  '^Z~'^X*^X^8t~ 


P^^X 

*  1“  ' 
cc 


and  (  )'=d(  )/dT,  (  )«=d(  )'/dT,  T=«t;  Al^^lt-ltc 

Tjj:  the  tension  of  the  tether 

Wj,:  the  orbital  angular  velocity 

nig:  the  mass  of  the  subsatellite 

m^:  the  mass  of  the  tether 

1^^:  the  command  length  of  the  tether 
mg^:  the  total  mass  of  the  subsatellite  and  the  tether 
h^;  coordinate  of  the  length  of  the  boom  in  the  shell 
reference  system 

Jg,  Jy  and  Jg  are  principal  moments  of  inertia  of  the 
shell 


If  we  let 


Af=3  (ic+p  )  +■ 


^  tc^  c^st 
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Q  _  sv* ^ tv 

®  2  T* 

(*)  cJy 


^  _  ^sx*^tx 
^  ec^  c^st 


then  the  equations  of  in-plane  notion  become: 


ky^a''+Q''*3a.-*-3d-k^e*c-D^ 


e!'t*2  (ic+P  jf)  Q^*2ka'-3ke 


The  system  equations  of  motion  can  be  written  in  matrix 

iC^=Ax+5u+CC 


where 


[0,a,€,0^,  a',e^]  ^ 


and 


o' 

0  I' 

,  B  = 

f°] 

B,  = 

0 

/ 

,1. 

It  is  obvious  that  if  we  consider: 

U=l^*Uc  >  ' 


(5.33) 

(5.34) 

(5.35) 

form  as: 

(5.36) 
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and  if  u_  can  be  found  such  that 


Bx^*cn=o 


then  the  system  equation  (5.36)  becomes 

,  (5.37) 


and  an  all  state  feedback  control  law  can  be  constructed: 

(5.38) 

For  the  system  parameters:  mj,=10,000  kg,  m3=500  kg,  111^=8. 35, 
l^g=l,000  m,  and  80  meter  boom,  an  optimal  control  u^  which 
minimizes  the  performance  index 

«• 

J-  f  (xI{2X-*-U^RUa)  dt  (5.39) 


can  be  found  by  carefully  selecting  Q  and  R. 
here  state  variable 

Q:  positive  semi-definite  state  penalty  matrix 

R:  positive  definite  control  penalty  matrix 

Some  typical  simulation  results  demonstrate  that  the  transient 
response  with  the  optimal  gains  obtained  by  suitable  Q  and  R  has 
low  overshoots,  short  settling  times  and  rising  times. 

But  the  optimal  control  based  only  on  the  use  of  tether  tension 
can  not  realize  Bji^+CC®0,  since  only  one  control  exists.  Thus,  the 
best  way  to  proceed  is  to  form  an  optimal  compensator: 

Lf^=-(B^B)-^B^CI2  (5.40) 

Nevertheless,  the  disturbance  D0  and  can  not  be  totally 
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compensated  in  this  way.  Fortunately,  the  effects  of  D9  and  are 
not  significant  in  lower  orbit.  On  the  other  hand,  the  in-plane 
motion  of  the  shell  and  the  tether  are  coupled  with  the 
longitudinal  motion  of  the  tether.  The  system  still  can  be 
controlled  to  a  satisfactory  degree  in  lower  orbit,  which  can  be 
demonstrated  by  some  numerical  simulation  results  (Fig.  5.7).  If 
the  altitude  level  is  increased  to  approximately  8600  km,  it  is 
difficult  to  maintain  the  satisfactory  pointing  accuracy  of  the 
shell  reflector  because  of  the  significant  effects  of  the 
uncompensated  disturbances.  Therefore,  in  order  to  still  control 
the  system  in  higher  orbit  with  the  presence  of  the  solar  radiation 
disturbances,  a  hybrid  control  system,  i.e.  introducing  some  other 
kind  of  actuator  on  the  shell  or  on  the  subsatellite,  is  necessary. 

5.4  inuierical  Siaulatioiis 

The  numerical  simulations  are  conducted  based  on  the  previous 
assvimptions.  Since  the  solar  incidence  vector  expressed  in  the 
orbital  reference  system  is  a  function  of  the  orbital  period  and 
the  motions  of  the  tethered  system  are  varied  periodically  for  the 
most  part,  the  disturbance  torques  induced  by  the  solar  radiation 
pressure  are  periodic  functions;  and  the  period  is  close  to  the 
orbital  period,  which  is  demonstrated  by  the  numerical  simulation 
results  (see  Figs.  5.3,  5.4,  5.5  and  5.6).  Comparing  Figs.  5.3, 
5.4,  5.5  and  5.6,  we  can  easily  see  that  the  magnitude  of  in-plane 
torques  are  much  larger  than  those  of  out-of-plane  torques,  and  the 
maximum  amplitude  of  the  in-plane  torques  in  the  completely 
reflecting  case  is  greater  than  that  in  the  completely  absorbing 
case.  In  the  completely  absorbing  case,  the  torque  acting  on  the 
subsatellite  about  the  axis  is  larger  than  the  torque  acting  on 
the  shell  about  the  axis,  but  vice  versa  in  the  completely 
reflecting  case.  Therefore,  if  the  solar  radiation  disturbances  are 
considered,  the  major  disturbances  are  the  in-plane  torques  acting 
on  the  shell  and  the  subsatellite  because  of  the  large  illuminated 


area  on  the  shell  and  the  long  moment  arm  from  the  center  of  mass 
of  the  whole  system  to  the  subsatellite. 

From  Fig.  5.7,  it  is  noticed  that  the  effects  of  the  solar 
radiation  disturbances  are  small  in  the  low  earth  orbit  (L.E.O), 
since  the  steady  state  performance  of  the  tethered  system 
controlled  by  the  tether  tension  with  the  solar  radiation 
disturbances  is  almost  the  same  as  that  without  the  solar  radiation 
disturbances.  But  in  the  high  earth  orbit  (H.E.O) ,  the  higher  the 
altitude  level  is,  the  worse  the  performance  of  the  steady  state  of 
the  tether  system  (see  Fig.  5.8);  therefore,  the  effects  can  not  be 
neglected,  especially,  when  the  altitude  level  is  higher  than  about 
8600  km.  This  suggests  that  the  attitude  and  dynamic  control  law 
for  such  a  system  in  a  high  earth  orbit  or  geosynchronous  orbit 
(G.E.O)  should  involve  some  kind  of  active  (actuator)  control  on 
the  shell  or  subsatellite. 

5.5  Conclusions 

Several  conclusions  can  be  obtained  from  the  numerical  simulation 
results: 

1.  The  effects  of  the  solar  radiation  disturbance  in  lower  orbit 
can  be  neglected  but  they  are  more  significant  in  higher 
orbits. 

2.  The  tension  control  law  with  optimal  control  gains  which  are 
obtained  by  carefully  selecting  the  state  penalty  matrix  Q  and 
control  penalty  matrix  R  is  able  to  maintain  the  satisfactory 
pointing  accuracy  when  the  altitude  is  lower  than  about  8634 
km,  for  a  100  m  diameter  shell  reflector  connected  to  a  1  km 
tether  at  the  end  of  an  80  m  boom. 

3 .  The  maximum  magnitude  of  the  solar  radiation  disturbance 
torq’-e  which  disturbs  the  in-plane  motion  of  the  system  is 
larger  than  that  which  disturbs  the  out-of -plane  motion  of  the 
system. 

The  dominant  contribution  to  the  solar  disturbance  torque  is 
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induced  by  the  solar  radiation  pressure  on  the  shell  reflector 
or  the  subsatellite. 

In  the  case  of  a  completely  absorbing  surface,  i.e.  6=0,  the 
magnitude  of  the  disturbance  torque  is  smaller  than  that  in 
the  case  of  a  completely  reflecting  surface,  i.e.  6=1. 

In  the  case  of  a  completely  absorbing  surface,  the  solar 
radiation  disturbance  torque  contributed  by  the  subsatellite 
is  greater  than  that  contributed  by  the  reflector  shell;  and 
vice  versa  in  the  case  of  a  completely  reflecting  surface. 

In  order  to  control  such  a  system  to  a  satisfactory  degree, 
a  hybrid  control  law  based  on  both  tension  modulation  and  some 
kind  active  (actuator)  is  needed. 
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Tethered  Antenna/Reflector  System 
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Mbz;  Torque  on  Boom  about  Zr  Axis 

Mtz:  Torque  on  Tether  about  Zt  Axis  Fig.  5.4 
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In-piane  Solar  Radiation  Torques 

Completely  Reflecting  Surface 


Mfy;  Torque  on  Shell  about  Yr  Axis 
Msy:  Torque  on  Subsat.  about  Yt  Axis 
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Tether  ControhftespeRse-  IH.E.G! 

With  Solar  Disturbance 


Altitude:  8634.0  km 
R-1 


Altitude;  5421.8  km 
R-1 


Fig.  5.8 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  momentum  exchange  controller  which  utilizes  the  time  rate  of 
change  of  the  flexible  momentum  relative  to  the  rigid  body  motion 
as  a  part  of  a  feedback  control  law  for  maneuvering  and  vibration 
suppression  is  proposed  and  analyzed.  For  this  initial  application 
a  model  of  a  rigid  hub  (base)  with  cantilevered  flexible  appendages 
undergoing  a  single^axis  maneuver  is  examined.  The  lower  and  upper 
bounds  of  the  control  gains  for  system  stability  are  obtained  based 
on  Lyapunov  methods.  For  further  vibration  suppression  an 
additional  independent  flexible  control  system  activing  on  the 
flexible  appendages  can  be  included. 

This  momentum  exchange  control  technique  has  also  been 
incorporated  in  an  adaptive  feedback  control  concept  for  the 
retrieval  of  the  subsatellite  in  the  orbiting  tethered 
antenna/ ref  lector  system.  The  time  rate  of  change  of  that  part  of 
the  angular  momentum  attributed  to  the  tether  motion  from  the 
local  vertical  is  used  as  a  part  of  the  feedback  law  in  the  antenna 
attitude  control  system.  This  is  used  in  conjunction  with  an  open- 
loop  exponentially  decreasing  commanded  length  during  retrieval. 
Numerical  simulation  results  verify  the  effectiveness  of  this 
approach  for  the  suppression  of  both  the  antenna  and  tether  in¬ 
plane  motions  at  quick  retrieval  rates. 

The  optimal  linear  quadratic  Gaussian  (LQ6)  digital  control  of 
tne  orbiting  tethered  antenna/ ref lector  system  is  analyzed.  In 
this  high  order  system  (up  to  40  state  components)  model  the 
flexibility  of  both  the  antenna  and  the  tether  are  included.  With 
eight  point  actuators  optimally  positioned  together  with  tether 
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tension  modulation  it  is  seen  that  the  degree  of  controllability  is 
very  high.  A  method  of  measuring  tether  transverse  motions  is 
proposed  and  is  required  to  guarantee  system  observability.  For 
such  a  dynamic  system  with  plant  and  measurement  noises,  the  LQG 
technique  is  effective  for  designing  the  best  combinations  of 
controller  and  observer  pole  locations  for  the  control  and 
estimation  processes.  The  discrete  time  data  system  is  more 
practical  for  the  design  of  a  controller  based  on  an  on-board 
computer  system. 

An  analytical  formulation  of  the  modelling  of  the  solar 
radiation  disturbance  on  a  class  of  large  tethered 
antenna/reflector  types  of  orbiting  structures  has  been  completed. 
The  effects  of  this  disturbance  during  station  keeping  operations 
have  been  evaluated  based  on  numerical  simulation  of  the  equations 
of  motion.  For  the  large  orbiting  100  m.  diameter  reflector 
considered  and  connected  to  a  1  km.  long  tether  at  the  end  of  an  80 
m  boom,  it  is  seen  that  the  effects  of  solar  radiation  in  the  lower 
orbits  can  be  neglected.  A  tether  tension  modulation  control  law 
where  the  gains  are  based  on  the  liner  quadratic  regulator  theory 
is  able  to  maintain  satisfactory  pointing  accuracy  when  the 
altitude  is  lower  than  8634  km.  The  dominant  contribution  to  the 
solar  disturbance  torque  is  induced  by  the  solar  radiation  pressure 
acting  on  the  (shell)  reflector  or  on  the  subsatellite.  For  the 
higher  altitudes  a  combination  of  tether  tension  modulation  together 
with  active  (actuator)  control  will  be  required. 

It  is  suggested  that  further  research  needs  to  be  conducted  on 
the  application  of  the  new  momentum  exchange  and  compensation 
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control  technique  to  flexible  and  telerobotic  systems,  beginning 
with  a  simple  two-link  manipulator  system  in  two  dimensions, 
followed  by  an  extension  to  the  three  dimensional  model.  As  a 
further  application  of  this  concept,  the  use  of  a  movable  rigid 
boom  connected  to  the  Shuttle  (or  reflector)  and  through  which  the 
tether  would  be  deployed/ retrieved  should  be  studied.  In  order  to 
prove  the  feasibility  of  the  tethered  antenna/reflector 
subsatellite  system,  test  scale  models  for  in-orbit  experiments 
need  to  be  further  studied.  The  incorporation  of  the  momentum 
exchange  and  compensation  control  technique  with  the  test  scale 
model  should  be  emphasized. 

In  order  to  fully  implement  a  control  system  for  the  tethered 
antenna /ref lector ,  the  design  of  a  reduced  order  controller  is 
suggested.  Both  direct  and  indirect  methods  for  designing  a  low 
order  controller  which  still  achieve  satisfactory  performance  and 
robustness  should  be  considered. 

Further  investigation  of  the  effects  of  solar  radiation  on  the 
tethered  antenna /ref lector  system  should  concentrate  on  the 
related  effects  of  solar  heating  and  the  practical  implementation 
of  the  required  control /measurement  systems.  Measurement  of  in¬ 
plane  and  out-of -plane  angular  displacements  as  well  as  the 
relatively  small  levels  of  the  differential  tension  required 
present  challenges  to  the  current  state-of-the-art  capability. 

Finally,  questions  related  to  the  validation  of  the  defensive 
rapid  retargeting  results  for  the  reflector  system  as  a  three 
dimensional  maneuver  can  be  addressed  by  disigning  some  experiments 
using  the  ASTREX  facility  at  the  Phillips  Laboratory  as  a  3-D 
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experimental  test  bed.  These  experiments  would  be  based  on 
mathematical  studies  of  rapid  retargeting  slews  of  the  ASTREX 
configuration,  first  based  on  a  suitable  reduced  order  model  of  the 
flexible  ASTREX  and  an  implementation  of  the  application  of 
Pontryagin's  maximum  principle  based  on  the  quasilinearization 
method  for  solving  the  associated  two  point  boundary  value  problem. 
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